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I. 



INTRODUCTION 



Recent developments in modem sea warfare, such as the stealth technology and 
new ECM techniques, mandate the use of passive means for the detection of targets, 
along with the use of radar. Infrared (IR) systems meet the requirements of such a passive 
sensor, and therefore are becoming an increasingly important part of integrated combat 
systems. 

The main restriction imposed on current IR sensors, such as Infrared Search and 
Track (ERST) systems or Forward Looking Infrared (FLER) systems, is that they have to 
operate in a heavily cluttered sea background. In the frequencies used by radar a detection 
is achieved when the signal to noise ratio is greater than one. Unlike the microwave 
wavelengths used in radar, in the BR wavelengths the target detection is more frequently 
based on the signal to background ratio, instead of the signal to noise ratio, because in the 
ER region of the electromagnetic spectrum the background is very intense. FLER systems 
particularly are contrast limited in performance rather than noise limited. 

Many techniques have been developed to deal with the problem of background 
suppression in IR sensors. These techniques, which include spatial filtering and advanced 
signal processing, improve the detection probability, but do not provide a satisfactory 
solution. Polarization filtering is a relatively new approach to the problem of improving 
the target contrast against sea background. 
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A number of studies, [Ref. 1,2, 3,4, 5,6, 7] in the past seven years by the Naval 
Postgraduate School, suggest that polarization effects in the observed IR radiation by the 
sea background could be exploited to improve target contrast. 

Atmospheric propagation prediction codes, such as LOWTRAN and MODTRAN, 
play an important role in the development of IR sensors. Thus far none of these codes 
included polarization effects. Recently, a research grade upgrade to MODTRAN 
developed by C.R. Zeisse at SPAWAR SYSTEMS CENTER [Ref. 8], has allowed the 
polarized case. The new code is the polarization version of SEARAD. It is capable of 
predicting the values of the polarized components of the observed sea IR radiance and the 
degree of polarization of the apparent radiance. 

This thesis uses experimental data available in the Electro-Optical Propagation 
Assessment in the Coastal Environment (EOPACE) data base, to perform a direct 
comparison to the SEARAD predictions. The specific data employed were taken by the 
Naval Postgraduate School during a multinational experiment in San Diego Bay on April 
1996. The data include IR images of the Research Vessel (R/V) POINT SUR, recorded 
using an AGA 780 camera. For a given period of time, the ship’s image against the sea 
background, was recorded with a polarizer filter in the horizontal direction, in the vertical 
direction, and removed. In that way, three sets of images were obtained: vertically 
polarized images, horizontally polarized images and unpolarized images. Meteorological 
and geographical data were recorded on board the R/V POINT SUR by the NPS 
Boundary Layer Meteorology Group and the NPS Center for Infrared Technology 
(NACIT). 
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The objective of this thesis is to correlate, and to analyze these data from the 
EOPACE database, and compare the results with those predicted by the polarized version 
of SEARAD, to validate the new code. 

This thesis begins in Chapter II with a brief discussion of ER fundamentals. In 
Chapter IH the polarization theory used by the SEARAD code is provided. Chapter IV 
describes the experimental data acquisition and analysis. Chapters V and VI present the 
results, the conclusions, and the recommendations, respectively. 
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II. INFRARED FUNDAMENTALS 



The analysis developed in this study deals with concepts related to infrared 
radiation. To better understand these concepts, this chapter presents the basics of infrared 
theory and terminology. 

A. THE NATURE OF INFRARED RADIATION 

The prefix “infra” has a Latin origin, meaning below or beneath. So the word 
infrared refers to the region of the electromagnetic spectrum below the red end of the 
visible light. Infrared radiation was discovered in 1800 by the English astronomer Sir 
William Herschel. All objects radiate infrared radiation, unless they have a temperature of 
0 Kelvins. IR radiation obeys Maxwell’s equations in its propagation. Furthermore, the 
same laws of diffraction, reflection, and refraction that apply to visible light, radio waves, 
and microwaves, also apply to the infrared radiation. 

The position of the infrared spectrum is between visible light and microwaves as 
shown in Figure 2.1. The wavelengths of infrared radiation (IR) vary from 0.7 pm to 1000 
pm. The infrared spectrum is further subdivided into IR bands, based on criteria such as 
the atmospheric transmission, the method of detection and the form of radiation. These 
bands are: 

• The near infrared (NTR) from 0.75 pm to 3.0 pm. 

• The mid wave infrared (MWIR) or Middle infrared from 3.0 pm to 6.0 pm. 

• The long wave infrared (LWIR) or Far infrared from 6.0 pm to 15.0 pm. 

• The extreme infrared (XIR) from 15.0 pm to 1000 pm. 
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Figure 2.1 Electromagnetic and Infrared Spectrum [Ref. 9] 

Reflected Solar radiation is the dominant component of infrared energy in the NIR 
band. On the other hand infrared radiation in the MWIR and LWIR bands, is 
predominantly due to object emission. The latter ER radiation is also called “thermal 
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radiation” because it is temperature dependent. Most practical applications of infrared do 
not include the XIR band. 

B. INFRARED TERMINOLOGY AND RADIOMETRIC UNITS 

The infrared community, like many other branches of engineering and science, 
uses its own terminology and units to describe the phenomena associated with the IR 
radiation. 

The basic units used are summarized in Table 2.1. Although the area units are 
frequently expressed in cm , in this thesis they will be expressed in m . The reason for 
this is that the polarized model of SEARAD uses meters squared as area units. 



Area (A) 


Square centimeters (cm 2 ) or Square meters 
(m 2 ) 


Wavelength (k) 


Micrometers (pm) 


Temperature (T) 


Degrees Kelvin (K) 


Plane angle 


Radian (rad) 


Solid angle (O) 


Steradian (Sr) 



Table 2.1 Basic Units Used in IR Theory 



Infrared theory uses specific radiometric quantities to describe the IR phenomena. 
These quantities are the following: 
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• Radiant Flux Density - Exitance (M): It is the radiant flux leaving an area of a 
surface divided by that surface 

• Irradiance (E): It is the radiant power incident on an area of a surface divided 
by that area 

• Radiant Intensity (I): It is the radiant power leaving a point source per unit 
solid angle. 

• Radiance (N): It is the radiant power leaving or arriving at a surface at a point 
in a given direction per solid angle and per unit area normal to that direction. 

The units used for the radiometric quantities are laid out in Table 2.2. As 
mentioned above, square meters are chosen as the area units, for compatibility with 
SEARAD. 

Often there is a need to refer to the radiometric quantities over a specific 
bandwidth and not over all wavelengths. When doing so the word “spectral” will precede 
the name of the radiometric quantity, and the corresponding units will be divided by the 
pm which is the unit of wavelength. 



Radiometric Quantity 


Units 


M 


Watts/m 2 


E 


Watts/m 2 


I 


Watts/Sr 


N 


Watts/m 2 Sr 



Table 2.2 Radiometric Units 
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C. RADIATION LAWS 

1. Blackbody Radiation and KirchhofFs Law 

The maximum power which can be radiated by an object at a given temperature is 
called blackbody radiation. A blackbody is an ideal object, which is both a perfect emitter 
and perfect absorber of radiation. An object which emits or absorbs less radiation than the 
blackbody, but with the same spectral distribution, is called a graybody. 

The emissivity s of an object is defined as the ratio of the exitance M e from the 
object to the exitance Mo from the blackbody at the same temperature: 



The absorptance a of an object can be defined as the ratio of the radiant flux per 
unit area being absorbed by the surface of the object to the radiant flux per unit area 
incident on the surface. 



Kirchhoff’s Law states that, if a graybody is placed in an enclosure in thermal 
equilibrium, then the spectral emissivity is equal to the spectral absorptance: 



£= Me/ M 0 



( 2 . 1 ) 



a= E a /Ej 



( 2 . 2 ) 



e(X)=a(X) 



( 2 . 3 ) 
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Equation 2.3 implies that the spectral radiant flux per unit area being emitted by 
the object is equal to the spectral radiant flux per unit area being absorbed by the object in 
a state of thermal equilibrium. This equation can be applied more widely to bodies not in 
thermal equilibrium. 

By the law of conservation of energy, the incident power P; upon an object, is 
equal to the sum of the transmitted power P T , the reflected power P p> and the absorbed 
power P a : 



Dividing both sides of the Equation 2.4 by Pi, and by defining x = P T /Pi as the 
transmittance, and p = Pp/Pi as the reflectance we get: 



For a blackbody there is no reflection or transmission i.e. x = p = 0, soa = e = 1. 
Also for opaque graybodies since there is no transmission, x =0: 



Pi- P T + P p + P a 



(2.4) 



x + p + a = 1 



(2.5) 



a = s = 1 - p 



( 2 . 6 ) 
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2 . 



Planck’s Law 



The blackbody radiates or absorbs, over a wide range of wavelengths. Max Planck 
first postulated that this radiation is due to vibrating atoms that emit or absorb radiation in 
discrete quanta of energy, which are given by the following equation: 



E = h v 



(2.7) 



where: 

• E is the energy change in Joules 

• h is the Planck’s constant equal to 6.6256xl0' 34 Joules sec 

• v is the frequency in Hz 

Planck’s radiation law, can be deduced from Equation 2.7. It was initially 
developed as an empirical fit to experimental measurements, where h is the fitting 
parameter, and is given by the following equation: 



M X CT) 



2nc 2 h/A. 5 

exp(-^-) - 1 
y XkT 



( 2 . 8 ) 



where: 

• Mx(T) is the spectral exitance at a temperature T, in W/cm 2 -pm. 

• T is the temperature of the source in K. 
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• X. is the wavelength in pm. 

• h is the Planck’s constant equal to 6.6256xl0° 4 Joules sec. 

• k is the Boltzmann’s constant equal to 1.38054x10^ J/K. 

O 

• c is the speed of light in vacuum equal to 3x10 m/sec. 

Planck’s law, gives the spectral exitance for a blackbody. The corresponding 
spectral exitance for a graybody can be obtained by multiplying the blackbody spectral 
exitance by the greybody’s emittance. Figure 2.2 gives the spectral exitance for a 
blackbody at certain temperatures. It is observed in this figure that the higher the 
temperature of the object is, the smaller the wavelength is, at which the peak of spectral 
exitance occurs. This is expressed by Wien’s displacement law described next. 

3. Wien’s Displacement Law 

By differentiating Planck’s Law with respect to wavelength and then by equating 
the result to zero we get the maximum of the spectral exitance for a specific temperature. 
Wien’s displacement law is the solution of the differentiation, and it is expressed as 
follows: 

X max T=2897.8 pmK (2.9) 

where: 

• ?'-max is the wavelength in pm, for which maximum exitance occurs. 

• T is the temperature of the source in K. 
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Figure 2.2 Graphical Representation of Planck’s Law [Ref. 9] 



4. Stefan - Boltzmann Law 

To obtain a closed form expression for the total exitance from a blackbody, 
Planck’s law can be integrated over all wavelengths. The result is known as the Stefan- 
Boltzmann law and it is given by the following equation: 
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( 2 . 10 ) 




where a is the Stefan-Boltzmann constant and is equal to 5.6697 xlO 12 W/cm 2 K 4 . 

For a graybody, which has an emissivity which is constant over all wavelengths, 
the exitance becomes: 



A practical way of calculating the spectral exitance without having to solve the 
Planck’s law integral directly is through the use of the universal blackbody curve depicted 
in Figure 2.3. This curve is a graphical representation of Planck’s Law expressed as a 
function of the product XT and normalized with respect to the peak wavelength. The 
integral referred to the right axis is the integral of the universal curve from zero to each 
value of XT. The “in band” spectral blackbody exitance is then given by: 



where f is the fraction of the exitance contained in a certain bandwidth and takes values 
from zero to one. 



M(T)=saT 4 



( 2 . 11 ) 



M(X,T)=fcT 4 



( 2 . 12 ) 
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Figure 2.3 Universal Blackbody Curve [Ref. 9] 

5. Lambert’s Law 

A blackbody with a rough planar surface emits radiation diffusely. Such a surface 
is called a Lambertian surface,, and the emitted radiant intensity is given by Lambert’s 
law, expressed by the following equation: 



1(0) = 1(0) cos(0) 



(2.13) 
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where 0 is the angle between the normal to the surface and the viewing line of sight, and 
1(0) is the normal radiant intensity. 

From Lambert’s law it can be shown that the radiance N radiated from a 
blackbody Lambertian source into a hemisphere is independent of the viewing angle, and 
has the value[Ref. 9]: 



N = M/ji 



(2.14) 



D. ATMOSPHERIC TRANSMITTANCE 

The laws described in Section C give the IR radiometric quantities for vacuum 
conditions. However, in most practical applications, the atmosphere is the medium 
through which ER radiation propagates. The atmosphere reduces the radiance reaching the 
sensor from the source by two main mechanisms: absorption, and non-forward scattering. 

The atmospheric absorption is due to individual molecules or aerosol particles 
such as fog, dust, salt, and others. The transmitted radiation is absorbed whenever it 
contains photons with energies that match the allowed transitional energies of the 
molecules or the aerosol particles. 

In general scattering can be classified as either forward scattering or non-forward 
scattering. The non-forward scattering decreases the received radiance, while the forward 
scattering increases the received radiance. The scattering occurs when the incoming 
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spectral radiation has wavelength components that are in the order of the size of 
molecules or aerosol particles. Then the incoming radiation is scattered in all directions. 

The atmospheric transmittance x accounts for the effects of radiance reduction by 
the absorption and the non-forward scattering. The atmospheric transmittance takes 
values from zero to one and can be calculated by the Lambert-Beer law as follows [Ref. 
10 ]: 



x (k)= exp(-p(A.)R) 



(2.15) 



where: 

• p(A.) is the spectral attenuation coefficient, or extinction coefficient which is a 
function of the wavelength X. 

• R is the path length. 

The atmospheric transmittance for a certain bandwidth is calculated by integrating 
Lambert - Beer’s law over the specific wavelengths of the bandwidth. This is done by all 
the standard atmospheric transmission codes such as LOWTRAN and MODTRAN. 

The atmospheric transmittance for one kilometer path range, a midlatitude 
location, and a rural aerosol, was calculated using LOWTRAN7 and plotted in Figure 2.4 
[Ref. 11]. In this figure we see that there are two regions of high atmospheric 
transmittance of the MWIR and LWIR bandwidths. These regions are called 
“atmospheric windows” and are from 3.0 pm to 5.0 pm and from 8 pm to 14 pm. Most of 
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the modem military infrared systems, such as Infra Red Search and Track (IRST) and 
Forward Looking Infrared (FLIR) operate in these two atmospheric windows of 
propagation. 




Wavelength U/m) 



Figure 2.4 Atmospheric Transmittance [Ref. 11]. 
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III. POLARIZATION OF SEA RADIANCE 



The phenomenon of the polarization of sea radiance has been addressed 
extensively in the existing literature [Ref. 1,2,3,4,5,6,7,12,13]. This chapter presents 
briefly the theory related to the polarization of sea radiance. The presentation is based on 
a technical report by C.R. Zeisse [Ref. 8], which accompanies the polarized version of 
SEARAD model that is evaluated in this thesis. 

A. POLARIZED REFLECTED RADIANCE 

In our preliminary analysis the reflection of monochromatic radiation from a flat 
sea surface will be considered. The sea surface is assumed to be tilted with respect to the 
horizon. The geometry of the problem is presented in Figure 3.1. In this figure the origin 
coincides with the point of reflection, the X axis is pointing in the up wind direction, and 
the Z axis is pointing to the zenith. The subscripts s and r refer to the source and the 
receiver respectively. The angle cp is the azimuth relative to the upwind direction and the 
angle 0 is the zenith angle. The vector U r points to the receiver, and U s points to the 
source. The unit vector normal to the sea facet is omitted. 

The electric vector of the incident and the reflected radiation has two components: 
the s component which is perpendicular to the plane of incidence and the p component 
which is parallel to the plane of incidence. The plane of incidence is defined by the unit 
vector normal to the sea facet, and the direction vector of the incident wave. 



19 



z 




Figure 3.1 The Coordinate System, 0 S is the Zenith Angle of the Source, 0 r is the Zenith 
Angle of the Receiver, and to is the Angle of incidence and the Angle of Reflection 

[Ref. 8]. 



The Fresnel equations can then be used to derive the reflected components of the 
electric field E 2 which are: 



E 2 s = P s exp(-i8 s )E 0 sin(a) (3.1) 

E 2p = p p exp(-i5 p )E 0 cos(a) 
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where: 



• Eo is the magnitude of the electric field of the incident radiation. 

• r is the amplitude reflection coeficient and r 2 = p is the reflectance. 

• 8 is the phase of the plane wave. 

• a is the angle of arrival (the angle between the E vector and the plane of 
incidence). 

The subscript s refers to the component of the electric field perpendicular to the 
plane of incidence, and the p refers to the component of the electric field parallel to the 
plane of incidence, referred to the facet. 

The polarizer has an axis called “polarization axis” along which the electric (E) 
vector passes. The components of the E vector in all other directions are blocked by the 
polarizer. If the polarization axis is parallel to the plane of the horizon, the horizontal 
component of the reflected electric vector along Uh in Figure 3.2 passes through the 
polarizer. On the other, hand if the polarization axis is normal to the plane of the horizon, 
the vertical component of the reflected electric vector along U v in Figure 3.2 passes 
through the polarizer. From Equations 3.1 and 3.2 and by the geometry presented in 
Figure 3.2, the horizontal E 2 h and the vertical E 2 V components of the reflected electric 
field are: 



E 2v =E 2p cos[3 + E 2s sin(3 ( 3 . 3 ) 

E 2h = E 2s cosP + E 2p sinP (3-4) 



21 



The angle (3 is the angle between the plane of incidence and the plane of vision. 
The plane of vision is the one containing the unit vector U r pointing to the receiver and 
the unit vector along the Z axis. 




PUNE OF 
INCIDENCE 



PUNE OF 
VISION 



Figure 3.2 Reflected Components of the Electric Field [Ref. 8] 
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The reflected radiance N r can be deduced from Equations 3.3 and 3.4, by taking 
the squared modulus E*E and assuming natural light as the incident radiation. Then we 
have the following equations for the reflectance p and its polarized components p h and p v 
[Ref. 8] 1 : 



= = — (cos P) 2 + — (sin P) 2 

N 0 2 2 



p . =^ = ^- (c °s|3 )2 +^- ( si n |3 ): 

JN 0 z 1 



(3.5) 



(3.6) 



P = Ph+Pv = 



N 0 




(3.7) 



where p p =r p and p s =r s are the reflectances for the p and s components, and N 0 is the 
incident radiance. 

B. POLARIZED EMITTED RADIANCE 

When we have radiation incident on the sea surface, part of that radiation will be 
reflected and part of the radiation will be transmitted into the water. If we consider the 
water to be thick none of the energy can escape by transmission. The energy is then 
absorbed by the water causing its temperature to rise. According to Kirchoff’s law the 
emissivity of the sea will be equal to the absorptance of the sea. The emissivity is equal to 



lln Ref. 8 the reflectance p is denoted by R. and it is referred to as the “Fresnel reflection coefficient”. 
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N sea /Nbb where N sea is the sea emitted radiance and Nbb is the blackbody radiance for the 
same temperature. The emissivity e is equal to the absorptance a, which is equal to the 
fraction of the radiant energy that the water absorbs from the incident radiation since 
there is no transmission: s = a = 1 - p ( Equation 2.6). This Equation also holds for the 
polarized components i.e. s s = a s = 1 - p s , and e p = a p = 1 - p p , 

Finally with an analysis similar to the one of the previous section, the following 
equations can be obtained for the emitted vertical N v sea , and horizontal N h sea polarized sea 
radiance, as well as the unpolarized N sea sea radiance [Ref. 8] 2 : 



NT 5 * 3 c £ 1 

£ h = ^ = y(cosp) 2 +-^-(sinP) 2 =--p h (3.8) 

^ bb ^ 



£ v 



XT 563 c £ 

= — — = — (sinP) 2 H — —(cos P) 2 
N bb 2 2 




(3.9) 



e = £ h +e v 



N sea 




(3.10) 



The emissivities e s , e p , and reflectances p s , p p , for sea water at a wavelength of 10 
pm, are shown in Figure 3.3. 



2 In Ref. 8 the emissivity e is denoted by A, and it is referred to as “the Fresnel emission coefficient”. 
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Figure 3.3 Sea Water Emissivities and Reflectances for a Wavelength of 10 pm; s Refers 
to the Component Perpendicular to the Plane of Incidence, and p to the Component 
Parallel to the Plane of Incidence [Ref. 8]. 



C. TOTAL RECEIVED RADIANCE 

The expressions derived in the two previous sections for the polarized reflected 
and emitted radiance refer to a single wave slope of the sea surface. Integrating over all 
slopes, and applying a weight for each one of them, it is possible to calculate the radiance 
leaving the footprint of the sea surface either by reflection or by emission. 

Usually the receiver is located a substantial distance away from the footprint. 
Thus atmospheric attenuation and path radiance must also be taken into account, along 
with the reflected and emitted radiation. The path radiance N path is unpolarized. 



25 



Averaging over all angles of arrival and assuming a perfect polarizer, half of the path 
radiance will be added to the vertical component, and half to the horizontal component of 
the received radiance [Ref. 8]. 

Summarizing the above, the received sea radiance, with the polarizer in the 
horizontal direction Nh, in the vertical direction N v , and removed N is: 



N h = 



N„ 



p s cos 2 p + p sin 2 p\ /l-p s cos 2 (3-p sin 2 p\ I N path f 



+ N 



bb \ 



xf +- 



( 3 . 11 ) 



N = 



N„ 



p s sin 2 P + p cos 2 P\ /l-p s sin 2 P~p cos 2 P' 



+ N 



bb \ 



r J 



Tf + 



]SfP ath f 



( 3 . 12 ) 



N = N h + N v = 



N 0 Ps ^ Pp ^ +N bb ^l- Ps ^ Pp ) [xf + ^i 



( 3 . 13 ) 



where: 

• No is the incident radiance. 

• Nbb is the emitted blackbody radiance of the sea water. 

• p s is the reflectance in the s plane (plane perpendicular to the plane of 
incidence). 

• p p is the reflectance in the s plane (plane parallel to the plane of incidence). 

• p is the angle between the plane of incidence and the plane of vision. 
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• x is the “in band” transmittance along the path 3 . 

• f is the “in band” fraction of radiance 4 . 



The angular brackets in Equations 3.1 1, 3.12, 3.13 represent the following integral 
over the slope space: 




II 8^^-p(C,.C y .w)dC x dC y 

co<7t/2 U n 

U r =const 



JJ 

CO <71/2 

U r =const 



CO SCO 

cos 6 



n 



P(Cx.C y . w )dC x dC y 



(3.14) 



where: 

• co is the angle of incidence and the angle of reflection of the radiation (Figure 
3.1). 

• 0 n is the zenith angle of the sea facet normal. 

• (x is the facet slope in the upwind direction (Figure 3.1). 

• £ y is the facet slope in the cross wind direction (Figure 3.1). 

• p(£ x > Cy> W) is the probability that a wave facet will have a slope within ±di^/2 
of ( x , and ±d( y /2 of ( y , with a wind speed W. This probability is given by a 
model developed by Cox and Munk [Ref. 14, 15] 



3 t is referred to as the “spectral transmission” in Ref. 8 

4 f is referred to as the “relative spectral responsivity” in Ref. 8 
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D. POLARIZATION CHARACTERISTICS OF SEA RADIANCE 

The radiance N 0 incident on the sea surface is the sum of radiance coming from 
the sun and radiance coming from the sky. Equations 3.11, 3.12, and 3.13 can be 
rewritten in a way such that the contribution of the sun incident radiation and the sky 
incident radiation is more explicit. By denoting N sun the integrated reflected radiance due 
to the sun radiation, N sky the integrated reflected radiance due to the sky radiation, and 
using the h and v subscripts, for the horizontal and vertical components respectively 
Equations 3.11, 3.12, and 3.13 can be rewritten as follows: 

N„={Nr+Nf + Nr}tf + |^|f (3.15) 

N.={Nr+Nf +Nf}tf + |i^-jf (3.16) 

N = N h +N V ={N sun +N sky +N sea }tf +{N path )f ( 3 . 17 ) 



The degree of polarization P of the observed sea radiance is then defined as: 



- N N - N 

h n v _ iN h 

N ~N h + N 



A negative value of the polarization indicates vertically polarized radiation. The 
percent degree of polarization is simply the degree of polarization multiplied by 100 . 
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The problem of predicting the polarization of the observed sea radiance has been 
addressed in previous studies [Ref. 1,2]. In these studies it has been shown that the 
reflected sun and sky radiance are predominantly polarized in the horizontal direction. On 
the other hand the emitted sea radiance has been shown to be vertically polarized. The 
total degree of polarization of the received sea radiance is given by summing the 
horizontally polarized reflection terms with the vertically polarized emission term. 

Figure 3.4 shows the predicted degree of polarization for the MWIR and the 
LWTR bands, in the sun glint region. The polarization is plotted, for selected sun zenith 
angles, as a function of the observer azimuth angle relative to the sun azimuth. In the case 
of the MWIR band, the polarization appears horizontal (positive) inside the sun glint, and 
vertical (negative) outside of it, for most sun zenith angles. In the LWER band, the 
polarization is vertical both within and outside the sun glint region, for almost all the 
selected sun zenith angles. The difference in polarization between the two bands is due to 
fact that the LWIR band contains only 0.1% of the vertically polarized reflected solar 
radiation [Ref. 2]. The SWIR band, on the other hand, contains 6% of the solar radiation 
[Ref. 2]; that is why the horizontally polarized reflected sun radiance dominates over the 
vertically polarized emitted sea radiance, in the sun glint region. 
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Experimental data analysis [Ref. 2] confirms the predicted degree of polarization 
depicted in Figure 3.4. Concluding it can be said that in the LWER band that is examined 
in this thesis, the observed sea radiance is polarized vertically outside the sun glint region, 
due to the dominant sea surface emission term of the received radiance. 





MWIR IWIR 

Figure 3.4 Theoretically Predicted Degree of Polarization in the Sun Glint Region, for the 
MWER and LWIR, as a Function of Solar Azimuth and Zenith Angles [Ref. 1], 
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IV. DATA ACQUISITION AND ANALYSIS 



To examine whether the polarized version of the SEARAD model predicts 
accurately the polarization effects of sea radiance, a direct comparison with experimental 
data is needed. Data are available to qualified people in the EOPACE data base in the 
internet address http://sunspot.spawar.navy.mi1/543//eopace/eomain.html. The acquisition 
of the data was done by Professor A.W. Cooper, W.J. Lentz, LT Britton, Dr. P.L. Walker, 
the crew of R/V POINT SUR, and the NPS Boundary Layer Meteorology Group by 
Professor K.L. Davidson and P. Fredrichsen. I performed the data analysis. 

A. EXPERIMENT DESCRIPTION 

The EOPACE data analyzed in this thesis were collected during an experiment 
conducted in April 1996 in San Diego Bay. Infrared images in the LWIR band of the R/V 
POINT SUR were recorded for two days, 9 and 10 of April, using the AGA-780 thermal 
imager. Internal polarization filters were employed with' the camera to obtain three sets of 
images, unpolarized, horizontally and vertically polarized, for a certain period of time and 
position of the ship. The camera was positioned on the tip of Point Loma, inside NRAD 
building 15, with geographical coordinates 32° 39’ 55” N, 117° 14’ 31” E. It was 
mounted on a pan and tilt head, extended through the window. The height of the camera 
was not available from the EOPACE data base. To overcome that, the known position of 
the camera, which was 32° 39’ 55.02” N and 117° 14’ 31.02” E, was plotted on a 
nautical chart of the San Diego Bay and the elevation was found to be 90 feet (27.4 m) 
from mean sea level, using the isoheights on that chart. The ship assumed predetermined 
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positions, with bearings ranging from 135° to 200°, and ranges from 0.5 to 3.0 nautical 
miles as indicated in Figure 4.1. In each of the above positions the ship changed its 
heading from 280° to 100° turning clockwise. Full meteorological data, including ship 
GPS position, heading, and ship temperatures, were recorded on the R/V POINT SUR, 
for the two days of the measurements. 




Figure 4 1 Experiment Location Map Depicting the Positions Assumed by R/V POINT 

SUR. 



32 



B. INSTRUMENTATION 



1. The AGA 780 Thermovision Thermal Imaging System 

The AGA 780 is a dual band infrared scanner. It can operate in the MWIR band 
from 3 to 5.6 pm using an InSb detector, and in the LWIR band from 8 to 12 pm using a 
HgCdTe detector. The cooling of the detectors is accomplished by the use of liquid 

o 

nitrogen at 77 K. The AGA scanner unit uses two rotating prisms, one vertical and one 
horizontal, producing one interlaced frame consisting of four fields each having 70 active 
image lines. There were two types of input lens used with the camera in the experiment, 
7x7 and 3.5 x 3.5 field of view lenses. The geometrical resolution achieved by the 
AGA 780 is 1.1 mrad for the 1 x 7° lenses and 0.5 mrad for the 3.5° x 3.5° lenses [Ref. 
16] The aperture was set to f/1.8 for the duration of the experiment. The arrangement of 
electro-optical components of the camera is presented in Figure A1 in Appendix A. 

The camera’s measure of the received radiance, is called the Thermal Value [Ref. 
16: p. 10.1]. It is measured in arbitrary Isotherm Units. The Thermal Value is proportional 
to the incident photon radiation flux. On the other hand, the relation between the Thermal 
Value and the apparent temperature of an object is not linear. The latter relation can be 
obtained by the means of laboratory calibration curves that give the Thermal Value as a 
function of blackbody temperature as shown in Figure A.2 in Appendix A. 

2. Internal Polarization Filters 

The polarization filters used in the experiment were manufactured by Graseby 
Specac, and consisted of an aluminum grid on a round KRS-5 substrate. The 
specifications of the filters predicted 85% transmission efficiency for wavelengths from 1 
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to 12 pm. The effectiveness of polarization was given as 98%. The performance curves 
for the two filters are depicted in Figures A. 3 and A.4 in Appendix A. As we can see from 
the two Figures the performance of the two filters is nearly identical. 

The filters were positioned in the internal filter wheel mounting of the AGA, 
which has seven positions, selectable from the filter control knob, as shown in Figure A.l 
in Appendix A. The positioning of the two filters was done in such a manner that one 
filter would pass the horizontal incident radiation, and the second one the vertical 
incident radiation. Thus by changing the position of the knob periodically, three sets of 
images of the R/V POINT SUR were recorded: 

• horizontally polarized images, in position 4 

• vertically polarized images, in position 5 

• unpolarized images, in position 6, with no filter. 

For each of the above positions a calibration procedure was performed in the 
laboratory, (by E.C. Crittenden on 3-21-96), using the NOSC LW Head. He recorded the 
AGA thermal level of a blackbody for various temperatures. Then by a curve-fitting 
process, a calibration curve was obtained for each of the three positions, 4, 5, and 6, as 
shown in Figure 4.2. The curve fitting equation is: 

1 = +OS t 4 - 1 ) 

C*exp(B/T)- 1 



where: 
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I = Thermal value corresponding to temperature T (Isotherm Units) 
T = Temperature in K 



• OS = Offset 



The constants A, B, C depend on the filter position set on the camera 4, 5, or 6, as 
shown on Figure 4.2 

Figure 4.2 Calibration Curves for Filter Positions 4,5,and 6. 



NRAD (NOSC) Head, lw 

Internal Polarizers. 3/21/96 




CAVITY BLACK BODY TEMPERATURE 
+ Position 5 o Position 6 x Position 4 

3. RA POINT SUR 

The R/V POINT SUR is a ship 135 feet long, which is owned by the National 
Science Foundation and operated by the Moss Landing Marine Laboratories. In the San 
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Diego experiment on April 96 it was used as an infrared target. The ship’s equipment was 
used to record meteorological data, and GPS position, and a thermistor set installed by 
NPS measured spot temperatures of its hull. 

C. ORIGINAL DATA 

1. Recording Software 

The experiment data were recorded using a program called PTRWIN developed 
by the French company CEDIP. This software employed five hardware cards, the PTR 
card, the DSP card, the RANGEL card, the TOPAZE card, and the VGA card, installed in 
the docking station of an NEC laptop computer for the recording of the data. These 
acquisition boards allowed the interface with the AGA 780. Of these cards the one 
labeled PTR by the manufacturer allowed digitization of the analog video of the scanner 
into 12 bit integers which are referred to as the Digital Level. Besides being used as a 
recording software, PTRWIN is also able to perform some limited analysis of the data 
under the Windows operating system. 

2. Images 

As mentioned above, three sets of images were recorded for a specific position 
and heading of R/V POINT SUR: horizontally polarized, vertically polarized and, 
unpolarized images. For each polarization case there is a series of images recorded, 
usually 26. A sample set of these images is contained in Appendix B. The polarized 
frames are recorded together in one file, and the unpolarized ones are recorded in a 
separate file. Each file begins with a main header which contains general information, 
such as day and time data, scanner parameters etc. The main header is then followed by 
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the image frames. Each frame consists by a frame header and the four fields of an image 
as shown in Figure 4.3 
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Figure 4.3 Structure of the Image Files 



Each field is made up of 64 lines. A line has 125 pixels. There are four fields per 
frame, so the frame is an array of 4(64)=256 by 125 pixels. 

3, Geographical and Meteorological Data 

The geographical data were recorded using the ship’s GPS system located in the 
electronics laboratory. For the collection of the meteorological data, two meteorological 
towers were installed on the ship. One was positioned over the bridge, and the other on 
the bow. The measurements from these towers were averaged over periods of 30 seconds 
and recorded along with the GPS position and heading into two different sets of files, the 
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METOC 1 and METOC 2 format files. The exact structure of those files is given in 
Appendix C 

4. External Ship Temperatures 

To measure the external temperatures of the ship, 16 thermistors were installed on 
her hull by NACIT- NPS. The thermistors were connected to a computer and their outputs 
were recorded along with time information in special ’’skin” files, for the duration of the 
experiment. The temperatures were recorded in 20 second intervals. The positions of the 
thermistors are shown in Figure D1 in Appendix D. The format of the files, in which the 
data were recorded, is also given in Appendix D. 

D. DATA ANALYSIS 

1. Data Processing 

I used programs written in Interactive Data Language (IDL) to process the data 
and bring them into a format suitable for analysis. 

a) Images 

The original image files with a PTW extension had to be converted to a 
different format with a PTE extension. The reason for that conversion was that the 
original files were too large, and contained information unnecessary for the purpose of 
this study. Furthermore the polarization cases were separated into two different files, one 
for the horizontal and vertical polarization and one for the unpolarized case. The IDL 
program used to perform this conversion, called CONVERT.PRO, is based on an existing 
program BASEGEN.PRO written by M.C. Pontes [Ref. 17]. 
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Another problem that arose with the original images was the fact that each 
pixel of the image was recorded in Digital Level. A pixel in the image occupies two 
bytes, but only twelve bits were used, giving Digital Level numbers from 0 to 4096. The 
Digital Level (DL) is an arbitrary unit of measurement which is linearly proportional to 
the received radiance. The existing calibration curves for the AGA 780 relate the 
isotherm units (IU) to temperature. Therefore, there is a need to convert from DL to IU, 
so that we can extract apparent temperature and radiance from an image. 

The relation between the DL and IU is not given by the existing PTRWIN 
documentation. There is however a feature in the analysis part of the PTRWIN program 
that allows the conversion from DL to IU for a single pixel. This is done by using 
information about the thermal range and thermal level settings of the AGA 780, at the 
time of recording. The thermal range is selected by a knob on the camera, as shown in 
Figure 4.4, and corresponds to different scales depending on the variation of intensity of 
the received radiation. The thermal level is also selected by the user in order to display 
both the lowest and the highest temperatures of the object. Both of these settings were 
recorded automatically for each frame in the original data, through the card labeled 
RANGEL by the manufacturer, in the NEC docking station. 
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Figure 4.4 Thermal Range and Level Controls [Ref. 16] 



The desired relation between DL and IU can be obtained by plotting the IU 
against DL for a large number of pixels taken from the same frame. Then a line fitted to 
the points gives the conversion relation for that frame, as we can see in Figure 4.5. 

After processing frames having different thermal ranges and levels it was 
concluded that a linear dependence exists between the DL and IU. This linear dependence 
is expressed by the following relation: 
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IU = a (DL) + b 



(4.2) 



where a and b depend on the thermal range and level as shown in Table 4.1. This process 
can be applied for each frame set used. The values a and b are constant between frames 
that were recorded with the same thermal range and level settings. 




Digital Level 



Figure 4.5 Plot of Isotherm Units against Digital Level. 
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Thermal Range 


a 


b 


5 


0.001221 


Thermal Level-2.5 


10 


2x0.001221=0.002442 


Thermal Level-2(2.5) 
= Thermal Level-5 


20 


4x0.001221=0.004884 


Thermal Level-4(2.5) 
= Thermal Level- 10 



Table 4.1 Values of the a and b Coefficients 



b) Geographical and Meteorological Data 

The original METOC data were recorded in GMT time, while the images 
were recorded in PST time, so the geographical and meteorological had to be converted to 
PST time, to be able to correlate them with a specific frame. Furthermore not all the 
information in the original files was useful for the present study, so a new file was 
created, called GEOMET.PTE containing only the parameters needed. Those parameters 
are listed in Table 4.2. 

c) Ship Temperature Data 

As mentioned in the previous section the ship’s external temperatures 
were recorded through the use of thermistors, in special “skin” temperature files. The first 
column, out of seventeen, of the original skin data files represents the time in seconds that 
has elapsed from the initial recording time as shown in Appendix D. There is not, 
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however, an indication inside the temperature files, or the files that accompany them, as 
to whether this time is local Pacific Standard Time or GMT time. 



Parameter 


Units 


Time elapsed from 00:00 PST 


Sec 


Latitude of the ship 


Min 


Longitude of the ship 


Min 


Antenna height 


Meters 


Air temperature 


Degrees C° 


Relative Humidity 


Percent % 


Air Pressure 


Mbar 


Wind speed 


meters/second 


Wind Direction 


Degrees relative to North 


Temperature of the Sea Surface 


Degrees C° 



Table 4.2 Parameters used in the GEOMET.PTE file 



To clarify this discrepancy, and be able to correlate the temperatures with 
the image frames, a direct comparison with the meteorological data was performed. 
Thermistor number three had been positioned in the water surface underneath the stem of 
R/V POINT SUR as depicted in Figure D.l. The temperature of the sea surface was also 
recorded by the meteorological stations. The two temperatures were plotted together for 
the 9th of April assuming that the thermistor temperatures were recorded in Pacific 
Standard Time (PST). This plot is shown in Figure 4.6. 
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The variation of the two temperatures over time, on Figure 4.6 is very 
similar, leading us to the conclusion that the ships hull temperatures were indeed recorded 
in PST time. 




o 5 10 is 20 25 

pst Time- 



Figure 4.6 Comparison of the Meteorological Data Sea Temperature and Thermistor 3 
Sea Temperature, to Clarify in which Time Zone were the “Skin “ Data Recorded 

2. IDL Programs 

The analysis of the experimental data was performed by using primarily programs 
written in the Interactive Data Language (IDL). IDL is an array oriented language which 
provides many mathematical analysis tools and techniques for graphical display. 
Compared to other programming languages such as FORTRAN or C, it is significantly 
faster because its operators work on entire arrays. Moreover it provides means of 
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visualization of the data. Finally its input and output capabilities allow the manipulation 
of any data file format. 

The purpose of the DDL programs written in this study was to take the original 
images, correlate them with the existing geographical, meteorological, and temperature 
data; and produce proper input for use by the SEARAD model. In that way a direct 
comparison can be made between the experimental data and the theoretical predictions of 
the polarized version of SEARAD, about the polarization effects from the sea surface. 

Some DDL programs of particular interest are described below. The code for these 
programs is given in Appendix E 

a) Geography. pro 

The SEARAD model requires as an input the relative geographical 
position of the observation point and the point on the sea we are interested in, which is 
called the footprint. Geography.pro takes as input the GPS position, and the height of the 
camera, as well as the GPS position, and the height of the ship which are recorded in the 
METOC data. Then based on that input, it calculates the slant path distance from the 
camera to the ship, the true azimuth of the ship, and the zenith angle of the line of sight to 
the ship, as shown in Figure 4.7. 

b) Data.pro 

Data.pro uses the recorded time on each image frame to interpolate the 
meteorological data, the position data given by Geography.pro, and the ship temperature 
data. After the interpolation a new file is created that has additional information about the 



45 



image including the thermistor temperatures. This additional information is inserted in 
the form of a matrix at the end of each image matrix. The contents of the matrix is laid 
out in Table 4.3. 



camera 



up 



True North 




Figure 4.7 Position Data given by Geography.pro 



c) Tape5ship.pro 

Tape5ship.pro reads the geographical and meteorological data created by 
Data.pro that correspond to a specific image frame. Then it produces a Tape5 file, which 
is an input file for the SEARAD code, for the position where the ship was for a specific 
time. The reason we need this file is that we want to know what the transmittance and the 
path radiance were for that position of the ship. We can get that information when we run 
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SEAR.AD with the Tape5 file that Tape5ship.pro produces, as input. A sample of the Out 
file that SEARAD creates is given in Appendix F. 



1 


Frame time elapsed in sec 


17 


Temperature from thermistor 1 in C 


2 


Latitude of the ship in min 


18 


Temperature from thermistor 2 in C 


3 


Longitude of the ship in min 


19 


Temperature from thermistor 3 in C 


4 


Ship height in meters 


20 


Temperature from thermistor 4 in C 


5 


Air temperature in C 


21 


Temperature from thermistor 5 in C 


6 


Relative Humidity in % 


22 


Temperature from thermistor 6 in C 


7 


Air Pressure in mbar 


23 


Temperature from thermistor 7 in C 


8 


Wind speed in m/s 


24 


Temperature from thermistor 8 in C 


9 


Wind Direction relative to North 


25 


Temperature from thermistor 9 in C 


10 


Latitude of the camera in min 


26 


Temperature from thermistor 10 in C 


11 


Longitude of the camera in min 


27 


Temperature from thermistor 1 1 in C 


12 


Height of the camera 


28 


Temperature from thermistor 12 in C 


13 


Slant path range from the camera to the 
ship in m 


29 


Temperature from thermistor 13 in C 


14 


Zenith angle from the camera to the ship 
in min 


30 


Temperature from thermistor 14 in C 


15 


True Azimuth angle of the ship in 
minutes 


31 


Temperature from thermistor 15 in C 


16 


Stack temperature given by the average of 
thermistors 12 and 14 


32 


Temperature from thermistor 16 in C 



Table 4.3 Information Matrix generated by Data.pro 
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d) Input. pro 

Input. pro can be considered as the main analysis program in this study. It 
employs many other programs -functions to extract input information for SEARAD, from 
a specific set of images which consists of a horizontally polarized image, a vertically 
polarized image and an unpolarized image. 

It begins by using the following programs, written by M.C. Pontes [Ref. 
17] and modified for the needs of this thesis: 

• Hotspot. pro: Finds the point on the ship’s stack with the greatest temperature 
(hotspot). 

• Horizon. pro: Provides the equation of the horizon in an image. 

• Elements.pro: Generates an array that has the same dimensions as the image 
array, but the values in this array are ones where there is a sea pixel, and zeros 
everywhere else (ship and sky). This array is called a mask. 

The program continues by sampling the selected image every nine pixels 
horizontally and every nine pixels vertically. This action corresponds to sampling the 
image every half degree vertically and every one fourth of a degree horizontally, because 
the image is an array of 125 by 256 pixels and the field of view of the camera is 7 x 7 . 
So 9 pixels vertically are equivalent to 7°/125(9) =0.5°, and 9 pixels horizontally are 
equivalent to 7 /256(9) =0.25 ' 

Input. pro then retrieves information about the zenith and azimuth angles of 
the reference point on the ship, and calculates the zenith and azimuth angles at every 
point that is sampled, by determining the relative position of the sampling point and the 
reference point. This is better understood by looking at an example. Suppose that at the 
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reference point the zenith angle is 91 , the azimuth angle is 100 , and the data coordinates 
are 50 and 130. One sample point has data coordinates 60 and 240. Then at the sample 
point the zenith angle is 91 + (60-50)/9(0.5 )= 91.56 and the azimuth angle is 100 - 
(240-130)/9(0.25) = 96.95°. 

Not all the sampled points correspond to sea pixels. To eliminate the 
unwanted points we multiply them by the mask matrix, obtained by the Elements. pro 
program. The mask matrix has ones where we have sea pixels and zeros everywhere else. 
After the multiplication process we are left solely with sea pixels, which are the ones we 
are interested in. 

Finally, Input.pro provides information about the wind speed, the wind 
direction relative to the LOS at each point, sea surface temperature, and time of each 
frame. These data along with the zenith and azimuth angles of each point are the only 
information needed for SEARAD to predict the radiance from a specific point at sea. 

e) Isotherm.pro, Radiance. pro, Invrad.pro, Offset.pro, and 

Temp.pro 

The task of providing SEARAD with inputs is accomplished by Input.pro 
as described above. However to be able to perform a direct comparison between the 
output apparent temperatures of SEARAD and the experiment, we need to have an 
apparent temperature for each of the sampling points. 

The data as mentioned earlier in this chapter were recorded in units of 
Digital Level. Isotherm.pro converts the digital level to isotherm units, using the linear 
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relations between the two, which depend on the thermal range and thermal level of the 
camera. 

The conversion from isotherm units to temperature can be achieved 
through the use of the calibration curves for each of the polarization filters. These curves, 
however, correspond to zero range measurements, so to use them, we must also bring the 
experiment measurements to zero range. The recorded isotherm units for each pixel in the 
image have the following components: 

• The radiance due to the temperature of the body multiplied by the 
transmittance for the range from the object to the camera. 

• The reflected radiance from the object surroundings multiplied by the 
transmittance for the range from the body to the camera. 

• The path radiance. 

For the case of the stack (hotspot) we know the temperature from the 
ship’s thermistors, and SEARAD can calculate the transmittance as well as the path 
radiance from the camera to the ship having as an input the Tape5 file created by 
Tape5ship.pro. 

Radiance.pro calculates the radiance emitted by the hotspot using the 
Stefan-Boltzmann law: 



M(T) = oT 4 



(4.3) 



where: 
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• M(T):is the total radiant exitance from a black body surface 

• a : is the Stefan-Boltzmann constant = 5.6697 xlO' 12 Watts-cm" 2 K' 4 

• T: is the black body’s temperature 

Then the radiance is: 

M 

N = — f (4.4) 

71 

where: 

• N: is the “in band” Radiance 

• f: is the “in band” fraction of radiance in the 8 to 12 pm range. It is obtained 
from the universal blackbody curve [Ref. 9]. 

The calculated black body radiance for the hotspot is then multiplied by 
the transmittance, and the path radiance is added to the result giving us the total radiance 
that the camera receives at a specific ship-camera range. In the previous procedure the 
reflected radiance originating from the object surroundings is not taken into account, 
because there is not enough information to do so. Nonetheless its contribution is not 
considered very big, because the stack is much hotter than its surroundings, has an 
emissivity close to one, and moreover, this reflected radiance is multiplied by the 
transmittance and thus its weight is greatly reduced. 

Invrad.pro uses the total apparent radiance calculated by Radiance.pro, and 
finds an apparent temperature by employing the Stefan-Boltzmann law once more. This 
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temperature is the apparent temperature of the stack which includes the effects of 
atmospheric absorption and path radiance. 

The process of correcting the stack temperature can be expressed by the 
following equation: 



T , WC.,T/»+N’* 

A final ~ V r ) 

1(5 



(4.5) 



where: 

• T final : is the corrected temperature of the stack. 

• Tjnitiai" is the temperature of the stack recorded by the thermistors. 

• f: is the “in band” fraction of Radiance in the 8 to 12 pm range. 

• a : is the Stefan-Boltzmann constant = 5.6697 xlO' 12 Watts-cm’ 2 K' 4 . 

• t: is the “in band” atmospheric transmittance. 

• N path : is the “in band” path Radiance. 

Offset.pro employs the corrected apparent temperature along with the 
known isotherm unit (IU) value on the stack to calculate the offset OS by using Equation 
4.1 as shown below: 



OS = I (4.6) 

Cexp(B/T) -1 

Since the offset has been calculated, the calibration curves for the AGA 
780 can now be used. With the IU and the offset known for every sampling point. 
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Equation 4.1 can be solved to give us the apparent temperature of every point. The above 
procedure is performed by Temp.pro. There is however an approximation involved. The 
range of the different points at sea to the camera is either smaller or greater than the 
distance from the camera to the ship. This difference in range affects the transmittance 
and the path radiance and subsequently the offset, giving a different offset for every 
range. Nevertheless to calculate these different offsets we need to have a reference 
temperature at every distance, such as the stack temperature. Unfortunately we do not 
have these reference temperatures, so we use the offset calculated for the ship-camera 
range, making the necessary assumption that the offset does not change significantly for 
the ranges involved. 



f) Tape5.pro, Out.pro 

Tape5.pro has the task of gathering information from Input. pro and other 
programs to write formatted Tape5 files, ready to be input by SEAR AD. Every sampling 
point in an image corresponds to a different Tape5 file. 

After SEARAD is run for the different Tape5 files a number of output files 
called Out are generated. These files are read by Out.pro which also makes the 
comparison of the SEARAD and experimental received sea radiance and degrees of 
polarization. These results will be presented in the next chapter. 

3. Polarized Version of the SEARAD Code 

The polarized version of SEARAD, developed by C.R. Zeisse at SPAWAR- 
SYSCEN NRAD, is a modification of the standard atmospheric propagation code 
MODTRAN 2. Besides the classic results of MODTRAN, the polarized version of 
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SEARAD predicts the polarization components of the received sea radiance. The way this 
is done is explained in Chapter HI, and summarized briefly by C.R. Zeisse, as follows 
[Ref. 8:p.l]: 

The Fresnel coefficients are applied to the reflection of unpolarized 
monochromatic radiation by a sea water capillary wave facet. The result 
for a single tilted facet is integrated over all possible tilts and wave 
numbers to predict the polarization dependence of radiance received from 
the wind-ruffled ocean surface. 

The SEARAD model uses an input file called Tape5 to compute the received 
radiance from the sea surface. An example of such a file is given in Appendix F. The 
Tape5 contains all the necessary meteorological and geographical information to calculate 
the sea radiance. Each line in this file is called a “card”. The description of the 
significance of each number in the cards is explained in the LOWTRAN 7 User 
Instructions [Ref. 18]. The differences of SEARAD from these instructions are contained 
in the “Notes” file which accompanies the polarized SEARAD model. 

The methodology followed for creating the Tape5 files, using DDL programs, was 
described above. Generally the following parameters were used: 

• Midlatitude Summer atmospheric model. 

• Slant path to space from the camera height to the sea. 

• Program execution in thermal radiance mode. 

• Navy maritime aerosol model. 

• Air mass parameter 3. 

• Wind speed and relative wind direction from METOC data 
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• Sea temperature from METOC data 

• Camera position, zenith and azimuth and time data, from the processed 
METOC data. 

After creating the necessary Tape5 files for all the sampled sea points in an image, 
a batch file was used to run the program sequentially. The approximate number of Tape5 
files per image was 160. 

After SEARAD is executed several output files are created. Those are the Tape6, 
Tape7, Tape8, Path, Sky, Sun, and Out files. For the purpose of the present study the Out 
file was chosen for the analysis, because it describes adequately the received sea radiance 
and its polarized components. An example of such a file is given in Appendix F. As can 
be seen in this file, the received sea radiance and the polarization are provided for zero 
and full range. The full range results include the path attenuation and the path radiance. 
Since the received experimental sea radiance was recorded at a certain range, the full 
range results of SEARAD were chosen for the comparison between theory and actual 
data. The results of this comparison are laid out in the next chapter. 
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V. 



RESULTS 



In this chapter all the experimental results will be compared directly, whenever 
this is possible, to the predictions of the theoretical model used by the polarized version 
of the SEARAD code. The comparison will include the apparent received sea radiance, 
the degree of polarization, and the azimuth dependence of the polarization. 

A. APPARENT RECEIVED RADIANCE 

The main concern when trying to compare apparent received radiance is to 
determine what the offset is in Equation 4.1, so that the conversion from IU to 
temperature can be accurate. Ideally we need a blackbody for every range and every 
polarization filter to obtain a point on the calibration curve, and thus acquire the offset. 
As mentioned in Chapter IV, we only have one such point, the ship’s stack. Therefore the 
offset is correct only for the range from the ship to the camera, and approximately correct 
for sea pixels in the same image at a distance around this range. 

Furthermore, in the case of images recorded with the use of polarizers we cannot 
use the calibration method analyzed in Chapter IV. This is because, in this case, the 
hotspot temperature corresponds to a polarized radiance which is a fraction of the total 
radiance. There is no method by which we can compute what this fraction is. The only 
way we could obtain the offset for the polarized cases, would be if we had a calibration 
curve at a specific range, for each of the polarizers used by the camera (vertical or 
horizontal). 
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Summarizing, due to calibration restrictions, the comparison between the 
observed sea radiance from the experiment and the calculated sea radiance from 
SEARAD was done only for the unpolarized case. A set of 34 images was used, in which 
the ship-camera distance varied from 0.5 to 1.5 nautical miles and the plots created are 
presented in Appendix FI. An example of such a plot is shown in Figure 5.1. The 
experimental and the SEARAD radiance are plotted against the zenith angle. The zenith 
angle was defined in Chapter IV in Figure 4.5. It is the angle between the Line Of Sight 
(LOS), which is the line from the camera to a specific point in the sea, and the line that 
points vertically from the position of the camera. Different zenith angles correspond to 
different slant ranges. The bigger the zenith angle the smaller the slant range. 

Figure 5.1, as well as all the other figures in this Appendix FI, was created by 
taking approximately 160 points from the unpojarized image and by running SEARAD as 
many times. Each time SEARAD was run the combination of azimuth angle, zenith 
angle, and relative sea wind was changed to correspond to the conditions of the 
experiment, and the target point that was chosen. All the other meteorological conditions 
were kept the same for each image. 

By examining Figure 5.1 and the other plots in Appendix FI, we observe that 
there is a good agreement between the experimental and the predicted sea radiance. In 
fact for this particular case, where the ship camera slant range is 1399 m, the average 
absolute error is 0.36 Wm' 2 sr'' (about 1.2% of the experimental radiance) and the 
standard deviation is 0.26 Wm' 2 sr '. 
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XSEARAD Radiance 
O Experimental Radiance 




Figure 5.1 Experimental and SEARAD Sea Radiance as a Function of the Zenith Angle. 

Different Values of Radiance for the Same Zenith Angle Correspond to Different 

Azimuth Angles. 

Figure 5.1 includes equal numbers of experimental and computed values at each 
zenith angle. We also observe that, for a specific zenith angle we have various values of 
the radiance. The different values correspond to different azimuth angles from the camera 
to the target point. The azimuth angle is also defined in Figure 4.3; it is the angle between 
the projection of the LOS to the earth and the true North. The variation in radiance is to 
be expected, because when the azimuth changes, the wind direction relative to the LOS, 
and the relative position of the sun both change, facts that change the total received 
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radiance [Ref. 1,14, 15]. It is important to the assessment of its accuracy, that SEARAD 
is actually sensitive to these kinds of changes. 

However, the experimental variation in radiance is bigger than the SEARAD 
variation in radiance for the same zenith angle. In Figure 5.2 the radiance is plotted as a 
function of azimuth relative to the solar azimuth. For clarity only two zenith angles have 
been selected. We see in this figure that the variation of the experimental radiance 
appears random, which may be due to data analysis errors. On the other hand the 
SEARAD predicted radiance generally decreases when we get away from the sun’s 
azimuth as one might have expected, because then the solar reflected component of the 
apparent sea radiance becomes less intense. From Figure 5.2 we also see that the variation 
of the SEARAD radiance is smoother in greater zenith angles (smaller distances), and 
bigger in small zenith angles (greater distances). 

In the other plots of Appendix FI the agreement between SEARAD and 
experimental results is not as close as in Figure 5.1. The reason for the increasing 
differences between the two could be attributed to the fact that an accurate camera 
elevation is not available. This factor affects the calculated zenith angles, and thus 
distances from the camera to the different points at sea. Different distances mean that we 
have different path radiance, and different path attenuation. It is guessed that the mistakes 
made in these two elements of the total sea radiance cancel each other in the case of 
Figure 5.1. That may be the reason the results in that case give such a good agreement 
between theory and experiment. 
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Finally another factor affecting the quality of the results is the accuracy of the 
offset calculation. It is known that the offset changes as the experimental conditions, such 
as the wind, air temperature, and time of day, change. That is why the offset was 
calculated on a frame by frame basis. Furthermore it is suggested from the analysis of the 
data that the offset may also be a function of distance. 



-X- SEARAD Radiance 
— e- Experimental Radiance 




Azimuth Relative to Solar Azimuth (deg) 



Figure 5.2 Experimental and SEARAD Radiance as a Function of the Azimuth Angle 
Relative to Solar Azimuth. The Radiance is Plotted for Two Different Zenith Angles, 

90.7, and 92.2 Degrees. 
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In Figure 5.3 the calculated offset is plotted against the slant path range from the 
camera to the ship. We see from this figure that the offset decreases exponentially as the 
distance increases. The apparent temperature T and the radiance increase as the offset 
decreases (Equation 4.5). So the apparent radiance increases with increasing distance 
(decreasing zenith angle). 

Note in Figure 5.1 and in Figures Gl.l through G1.5 that the SEARAD radiance 
varies with zenith angle which corresponds to a distance change, while the experimental 
radiance almost levels out at about 91 . These figures were plotted assuming the offset 
does not change in a particular frame. If based on the analysis of the previous paragraph 
the offset does change with distance then the plotted experimental radiance would be 
bigger in the small zenith angles and smaller in the bigger zenith angles. This would 
change the shape of the curve for the experimental radiance and probably make it look 
like the shape of the SEARAD radiance. Unfortunately, as mentioned above the offset is 
also a function of the experimental conditions; thus it was not possible to extract a 
reliable equation for the offset as a function of zenith angle, because the data points in 
Figure 5.3 were extracted from different images recorded at various times of the day. 
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Figure 5.3 Plot of the Calculated Offset as a Function of the Ship Camera Distance. Data 
Points were Extracted for Different Scenes and Experimental Conditions. 

B. DEGREE OF POLARIZATION 

The degree of polarization P is, as discussed in Chapter ID, the ratio of the 
difference of the horizontal and the vertical components of the radiance, divided by the 
total radiance. To get the degree of polarization per cent, P', we multiply the degree of 
polarization by 100: 



P 





(5.1) 



where: 

• N h is the horizontal component of the radiance 

• N v is the horizontal component of the radiance 
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• N is the total received radiance 

To determine the polarization from SEARAD, the Out file was used, an example 
of which is contained in Appendix D. The Out file explicitly determines what the vertical 
and horizontal components of the radiance are, as well as what the total radiance is. In the 
case of the experimental data we do not have the radiance recorded, but we do have the 
IU values, which are directly proportional to the radiance. By assuming that the set of 
vertically and horizontally polarized images are recorded simultaneously the degree of 
polarization per cent is calculated as follows: 



p=ioo-^ 113 v 



IU h +IU v 



(5.2) 



In Equation 5.2 it is assumed that the total unpolarized sea radiance is the sum of 
the horizontal and vertical components of the sea radiance [Ref. 8:p.l9]. It is mandatory 
to make that assumption, because the IU values in the unpolarized image have a different 
offset from the polarized images. Therefore the unpolarized IU can not be compared 
directly with the IU in the vertically and horizontally polarized images, which have equal 
offsets, as is shown in Figure 4.2. 

Taking approximately 160 points in each polarized image, that have different 
zenith and azimuth angles, and by running SEARAD for each of these points, we obtain 
plots like the one presented in Figure 5.4 where the degree of polarization per cent is 
plotted against the zenith angle. The different values of degree of polarization, observed 
in Figure 5.4, for the same zenith angle correspond to different azimuth angles. Negative 
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values of polarization correspond to vertically polarized radiation. The full range of plots 
for different ship-camera slant ranges is presented in Appendix F2. 

In Figure 5.4 the average absolute error is 0.17 % and the standard deviation is 
0.26 %. The results for the average absolute error and the standard deviation, in the cases 
included in Appendix F2, give the same order of agreement between SEARAD and 
experiment. It is therefore concluded that the polarized version of SEARAD is a reliable 
model for calculating the sea polarization. 

Note in Figure 5.4 that the degree of polarization is reduced at small zenith angles 
which correspond to greater distances. The degree of polarization decreases with 
increasing distance because the polarized reflected and emitted components of sea 
radiance are attenuated through the atmosphere, and unpolarized path radiance is 
accumulated in the apparent radiance causing the degree of polarization to be reduced 
[Ref. 8], 

Combining the results of this section with the results of Section A, the 
experimental received radiance for either the vertical or the horizontal polarization can be 
calculated [Ref. 1,2]: 



N h = N 



1 + P 



(5.3) 



where: 



N„ 



= N 



1-P 

2 



(5.4) 
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• Nh is the horizontal component of the radiance 

• N v is the vertical component of the radiance 

• N is the total received radiance 

Using the above equations, the problems of the polarized calibration can be 
overcome, provided that SEARAD is considered reliable. 



O Experimental Degree of Polarization (%) 
X SEARAD Degree of Polarization (%) 




90 90.5 91 91.5 92 92.5 93 93.5 

Zenith Angle (deg) 

Figure 5.4 Comparison of the SEARAD and Experimental Degrees of Polarization 
Different Values of Polarization for the Same Zenith Angle Correspond to Different 

Azimuth Angles. 
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C. AZIMUTH DEPENDENCE OF POLARIZATION 

The results presented in the previous sections indicate that there is an azimuth 
dependence of the degree of polarization of the total received sea radiance. Using the 
same assumptions and the same procedures as in the previous section, the experimental 
and the SEARAD values of the degree of polarization were plotted. This time instead of 
plotting the degree of polarization against the zenith angle, it was plotted against the 
azimuth angle relative to the sun. This angle is defined as the angle between the true 
azimuth of the LOS and the true azimuth of the sun. The angle calculation was based on 
the subsolar longitude and latitude given by SEARAD in the Out files (Appendix F), the 
camera’s geographical position, and the position of the sampled points of the sea surface. 
An example of such a plot is Figure 5.5 

It can be seen from Figure 5.5 that, for the chosen zenith angle, SEARAD predicts 
that the polarization will change slightly when going away from the sun direction, and it 
will become more negative i.e. more vertical. The fact that the received radiation is less 
vertically polarized in the region close to the sun has been observed in previous studies 
[Ref. 1,2]. In the region of sun glint the horizontally polarized solar reflected component 
of the total sea radiance becomes stronger as explained in Chapter IH. Thus the observed 
radiance is less vertically polarized when approaching the sun azimuth. 

The experimental degree of polarization, however does not seem to change as 
regularly as the one from SEARAD does. In fact it appears to be random. The sea 
footprint, i.e. the patch of sea surface that the camera sees, is 2.8 km or 0.82 nm“ for the 
image examined in the two previous sections. The sea footprints, for the rest of the 
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images are in the same order of magnitude. Based on the size of the footprint we may 
conjecture that this irregularity of the experimental polarization is attributable to the fact 
that the wind direction and strength are not correlated in the entire scene of the footprint. 
A change in the wind parameters changes the probability that a sea wave facet will occur 
within a specific slope, given by Cox and Munk [Ref.14, 15], and subsequently the final 
polarization of the observed radiance. (Equations 3.11, 3.12, 3.13, 3.14). 



—©—Experimental Degree of Polarization (%) 
X SEARAD Degree of Polarization (%) 




Azimuth Relative to Solar Azimuth (deg) 
Figure 5.5 Azimuth Dependence of Polarization 



In the analysis of the experiment, when calculations were made using SEARAD 
the conjectured variation of the wind parameters was not taken into account. The reason 
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for assuming the same wind conditions throughout the image is that the meteorological 
data were only recorded in the ship’s immediate area. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The primary research objective of this study was to determine whether the 
polarized version of the SEARAD model predicts accurately the apparent sea radiance, 
and the related degree of polarization, when compared with the EOPACE experimental 
data. 

The comparison of the unpolarized radiance calculated by SEARAD, with the 
experimentally recorded radiance showed that there is approximate agreement between 
the two, within 5%. More specifically the total average absolute error for all the selected 
images was 1.57 Wm‘ sr which is approximately 5 % of the experimental radiance. This 
error may be attributed, in part, to measurement uncertainties and calibration limitations. 
Overall the level of agreement between the SEARAD code and the experimental 
unpolarized radiance is considered satisfactory. 

The degree of polarization of the received sea radiance predicted by SEARAD 
agrees to a great extent with the value extracted from the experimental data. The total 
average magnitude of the difference between the two Pc^^-Pc was 

only 0.51 %. The results in this case were not affected by calibration limitations as in the 
radiance case. One possible cause of the error may be the fact that the selected 
experimental images were not recorded simultaneously, but with a delay of about five 
seconds one from the other. 

Summarizing, the analysis of the experimental data from the EOPACE database 
shows that the polarization version of SEARAD adequately predicts the received sea 
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radiance and its polarized components. Thus the polarization version of the SEARAD 
model is validated. 

The results of this thesis indicate that the observed sea radiance in the LWBR is 
polarized in the vertical direction, outside the sun glint. Therefore the results from 
previous studies by NPS, on the subject of the polarization of the sea radiance [Ref. 1, 2, 
3, 4, 5, 6, 7], indicating a suppression of the sea background through the use of polarizers 
are validated. 

The plots of the experimental and the SEARAD degrees of polarization suggest 
that the degree of polarization is reduced when the distance between the sea footprint and 
the sensor increases. Thus the contrast improvement with the method of polarization 
filtering of the sea radiance is maximized in small distances. 

In conclusion, the polarization version of SEARAD is considered a reliable model 
for the further study of the polarization effects of the sea radiance. This model could be 
used to identify the conditions under which polarization filtering gives the maximum 
background suppression and contrast improvement. A new experiment may then be 
designed to fit those particular conditions when that is possible. SEARAD could also be 
used in the development of a new generation of IR systems exploiting the polarization 
effects on the sea environment. 
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APPENDIX A 



AGA 780 THERMOVISION AND POLARIZATION FILTERS PERFORMANCE 

CURVES 

Figure A1 depicts the arrangement of electro-optical components inside the AGA 
780 [Ref. 16]. Figure A2 provides typical calibration curves for the camera, using various 
apertures, in Isotherm Units versus Temperature. 

Figures A3 and A4 present the performance curves of the two polarization filters 
used in the EOPACE April 96 experiment. In each of these Figures there are three curves. 
Curve one provides the transmittance of the filter with the grid normal to the incident 
electric field. Curve two shows the transmittance of the filter with the grid parallel to the 
incident electric field. Finally curve three gives the transmittance when the two filters are 
placed perpendicular to each other (i.e., crossed). The horizontal axis shows the 
wavelength and the vertical axis the transmittance. These curves were provided by the 
manufacturer. 
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Figure A.l Arrangement of electro-optical components inside AGA 780. [Ref. 16] 
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ISOTHERM UNITS 




Figure A. 2 Calibration curves for the AGA Thermovision 780 in the LWIR. [Ref. 16] 
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Figure A.3 Performance Curves for the First Polarization Filter 
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Figure A.4 Performance Curves for the Second Polarization Filter 
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APPENDIX B 



IMAGES 

Three images are presented in this Appendix. In the first image the received 
radiation is horizontally polarized, in the second one the received polarization is vertical, 
and in the third one the received radiation is unpolarized. Note the improved contrast of 
ship over sea in Figure B.l. Figure B.l shows much higher contrast than Figure B.2 or 
Figure B.3. because the background sea radiance is vertically polarized and by filtering 
out the vertical component of the radiance the contrast is improved. 
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Figure B.l AGA 780 Image with the Polarizer in the Horizontal Position. The Hotspot 
Appears as a White Spot in the Image. Progression Through the Color Spectrum Indicates 

Increasing Temperature. 
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Figure B.2 AGA 780 Image with the Polarizer in the Vertical Position. The Hotspot 
Appears as a White Spot in the Image. Progression Through the Color Spectrum Indicates 

Increasing Temperature. 
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Figure B.3 AGA 780 Image with the Polarizer Removed. The Hotspot Appears as a 
White Spot in the Image. Progression Through the Color Spectrum Indicates Increasing 

Temperature. 
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APPENDIX C 



METEOROLOGICAL FILES 



This Appendix presents the format description of the meteorological - 
geographical files used in the analysis of this study. The files were recorded in two 
formats, the METOCOl, and the METOC02 format. The storage format description 
included in this Appendix was provided by the NPS meteorological department. 
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STORAGE FORMAT FOR METOCOl (DASC) DATA 



The data is stored in a set of files with the names "eopace.*" and "30sec.*".The extension 
for the file names is the Julian date the file was started on. A new pair of files is started every day 
at 00:00 GMT. The time data in the files are all GMT. The eopace files have the data averaged 
over 10 minutes. The 30sec files have the GPS data averaged over 30 seconds and the 
last set of data from the Campbell. In the 30sec file, the Campbell data are not averaged. The 
data format for both files is the same. Each file has a header which tells when the file was started. 
Each file also has a footer that tells when the file was finished. Some files may also have 
additional start and stop times interspersed with the data, depending on whether data logging was 
turned off or the program was stopped before the end of the day. The data are arranged in fields 
separated by spaces across a line of text. A carriage return character and a line feed character are 
at the end of each line. Each new reading is on a new line. The fields are in the order below. The 
numbers before the fields are not in the data. They serve to make identifying the fields easier. 

01 number of readings in this average, nn 

02 year, yyyy 

03 Julian date, ddd 

04 hours & minutes, hhmm 

05 seconds, ss 

06 Relative Wind Speed, m/s, mmm 

07 Relative Wind Direction, degrees, ddd 

08 Ship Speed, knots, at present always 0 

09 Ship Direction, degrees, at present always 090 

10 True Wind Speed, m/s, mmm 

1 1 True Wind Direction, degrees, ddd 

12 T air, degrees C, tt.t 

13 RH, per cent, pp 

14 Pressure, millibars, pppp.p 

15 Sea Surface Temperature, degrees C, tt.t 

16 GPS time, hhmmss 

17 GPS Latitude, ddmm.m 

18 GPS North/South Indicator, either ‘N’ or 'S’ 

19 GPS Longitude, dddmiam 

20 GPS East/West Indicator, either 'E or ’W 

21 GPS Speed Over Ground, knots, kk.k 

22 GPS Course Over Ground, degrees, ddd.d 

23 GPS Antenna Height, meters, mm.m 

Here is a sample of the data: 

30 1996 99 1103 28 5 009 0 090 3 323 12.3 91 1013.9 

20.3 110838 3306.2 N 11806.0 W 0.0 90.0 189.6 
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STORAGE FORMAT FOR METOC02 DATA 



The METOC02 system had two additional temperature sensors and two additional RH sensors. The 
data are arranged slightly differently than the METOCOl system. The data is stored in a set of files with the 
names "eopace.*" and "30sec.*". The extension for the file names is the Julian date the file was started on. A 
new pair of files is started every day at 00:00 GMT. The time data in the files are all GMT. The eopace files 
have the data averaged over 1 0 minutes. 

The 30sec files have the GPS data averaged over 30 seconds and the last set of data from the Campbell. In 
the 30sec file, the Campbell data are not averaged. The data format for both files is the same. Each file has a 
header which tells when the file was started. Each file also has a footer that tells when the file was finished. 
Some files may also have additional start and stop times interspersed with the data, depending on whether 
data logging was turned off or the program was stopped before the end of the day. The data are arranged in 
fields separated by spaces across a line of text A carriage return character and a line feed character are at the 
end of each line. Each new reading is on a new line. The fields are in the order below. The numbers before 
the fields are not in the data. They serve to make identifying the fields easier. 

0 1 number of readings in this average, nn 

02 year, yyyy 

03 Julian date, ddd 

04 hours & minutes, hhmm 

05 seconds, ss 

06 Relative Wind Speed, m/s, mmm 

07 Relative Wind Direction, degrees, ddd 

08 Ship Speed, knots, at present always 0 

09 Ship Direction, degrees, at present always 090 

10 True Wind Speed, m/s, mmm 

1 1 True Wind Direction, degrees, ddd 

12 T air, degrees C, tLt 
13RH,percent,pp 

14 Pressure, millibars, pppp.p 

15 Sea Surface Temperature, degrees C, tt.t 

16 T air A, degrees C, tt.t 

17 RH A, per cent, pp 

1 8 T air B, degrees C, tt.t 
19RHB, per cent, pp 

20 delimiter between Campbell data and GPS data, "***" 

2 1 GPS time, hhmmss 

22 GPS Latitude, ddmnun 

23 GPS North/South Indicator, either 'N’ or 'S' 

24 GPS Longitude, dddmnun 

25 GPS East/West Indicator, either TS' or 'W 

26 GPS Speed Over Ground, knots, kk.k 

27 GPS Course Over Ground, degrees, ddd.d 

28 GPS Antenna Height, meters, mm.m 

Here is a sample of the data : 

18 1996 99 0319 28 2 090 0 090 3 263 12.2 96 1012.3 16.4 

12.2 96.9 12.7 95.7*** 031908 33223 N 11836.9 W 5.0 126.3 25.7 
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APPENDIX D 



THERMISTOR LOCATIONS AND TEMPERATURE FILES 

This Appendix presents the location of the thermistors placed on R/V POINT 
SUR, to measure the external “skin” temperatures of the ship. The temperature data from 
the 16 external thermistors were recorded by NPS - NACIT in special “skin“ temperature 
files along with time information. The description of the format of these files, and a 
sample of them, is also contained in this Appendix. 
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THERMISTOR LOCATION 




Flir 13 April 1, 1996 
PTSur 

Channel # Location 



Reference 25 0 C 
Aft Port 

Aft water probe 
Aft of Stack 
Aft On Stack 
Aft Boathouse 
Aft Starboard 
Port Boathouse 
Reference 25.0 C 
Port Below Stack 
Port Bow 
Starboard Stack 
Air 

Port Stack 
Bow Boathouse 
Starboard Boathouse 
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“SKIN” TEMPERATURE FILE FORMAT 



Data Format: Header - self explanatory Time in seconds, Channel 1 reference (25.001, 

channel2,3,4,5,6,7,8,Channel 9 reference (25.001), 10, 11,12,13,14,15,16 

aprill.l 

Mon,Day,Year= 4 1 1996 
Hour, Min, Sec= 20 8 6 

volts, istart,istop,nsam,iref 1.49939 1 16 100 1 

72486 25.001 12.863 14.185 12.883 14.996 14.756 12.581 13.813 25.001 13.034 13.182 
18.452 18.542 13.949 13.425 13.243 

72506 25.001 12.865 14.188 12.886 15.001 14.757 12.584 13.813 25.001 13.038 13.187 
18.461 18.546 13.951 13.423 13.246 

72526 25.001 12.868 14.140 12.888 15.006 14.756 12.588 13.817 25.001 13.039 13.189 

18.432 18.551 13.953 13.425 13.245 

72546 25.001 12.873 14.113 12.892 15.005 14.756 12.588 13.820 25.001 13.043 13.197 
18.445 18.557 13.958 13.424 13.243 

72566 25.001 12.877 14.118 12.896 15.003 14.755 12.590 13.822 25.001 13.045 13.201 

18.433 18.560 13.963 13.424 13.240 
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APPENDIX E 



IDL PROGRAMS 



This Appendix presents some of the most important DDL programs used in the 
analysis of the experimental data. 
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GEOGRAPHY.PRO 



.************** LTJG PANAGIOTIS KARA V AS HELLENIC NAVY***************** 

.********************** thesis research** ***************** ************ 

;*** VALIDATION EVALUATION OF THE POLARIZATION VERSION OF SEARAD****** 
;PROGRAM NAME :Geography.pro 

; Purpose: Calculates position of the ship relative to the camera 
function Geography, gpsship, gpscam 
; dimensions & initial values 
auxl = size(gpsship) 
aux2 = size(gpscam) 

; return if error 

if (auxl(2) ne aux2(2)) then return, input error’ 

; dimensions & definitions 
ncol = 3 
nlin = auxl (2) 
latship = gpsship(0,*) 
lonship = gpsship(l,*) 

hship = 7.3 estimated height above sea level of the ship 

latcam = gpscam(0,*) 

loncam = gpscam( 1 ,*) 

hcam = gpscam(2,*) 

pi = 3.1416 

fi = fltarr(l,nlin) 

; latitude variation = latidude distance in nmi 
dlat = latship-latcam 

latdist = dlat* 1852.0 ; convert nmi to m 

; longitude variation 
dlon = lonship-loncam 
; altitude variation 
dh = hcam- hship 
; midlatitude 

mlat = (latship+latcam)/2.0 
; transform midlatitude to radians 
mlat = mlat*pi/(60.0* 180.0) 

; longitude distance 

londist = dlon*cos(mlat)*1852.0 ; converted nmi to m 
; calculate fi angle 

for i = 0, nlin-1 do if (dlat(i) ne 0.0) then fi(i) =(-l)*atan(londist(i)/Iatdist(i))*(dlat(i) gt 0)+((- 
1)* atan(londist(i)/latdist(i))+pi)*(dlat(i) It 0) else fi(i) = (pi/2)*(dlon(i) le 0) + (1.5*pi)*(dlon(i) gt 
0) 

; calculate range R 
R =sqrt(dh A 2+latdist A 2+londist A 2) 

; calculate theta angle 
theta = pi-acos(dh/R) 

; return 

saida = [R,theta,fi] 
return, saida 
end 
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DATA.PRO 



.************** LTJG PANAGIOTIS KARAVAS HELLENIC NAVY***************** 
.********************** THESIS RFSF. A~R PH** ***************** ************ 

;*** VALIDATION EVALUATION OF THE POLARIZATION VERSION OF SEARAD****** 
;PROGRAM NAME :Data.pro 

;Purpose : Interpolate meteorological and geographical data into images creating a new image 
;array with additional information 
pro Data 

basepath =’c:\Thesis\data\basicV ; where is the BASEN.PTE 

filename = pickfile(/read, path = basepath , filter = ^.pte’, title=Select file to Analysis’) 

; Get Information 

info = leia_pte(filename, r info’,0) 

npol = info(7) 

ncol = info(8) 

nlin = info(9) 

nfra = lonarr(npol) 

for i = 0, npol - 1 do nfra(i) = info(10+i) 

; Size of elements (variables) in *.p te (in bytes) 



; offset for each result matrix (case) 
; 2 bytes per integer element 



; 4 bytes per float element / 100 columns 



resuoff = lonarr(npol) 
stime = nfra *2 *4 
slrange=nfra*2*2 
simag = nfra*ncol*nlin*2 
sresu = (nfra+l)*4*100 
; offsets 

timeoff = long(40+2*npol) ; offset for first time matrix 

lrangeoff=long(timeoff + total(stime)) 

imagoff = long(lrangeoff + total (slrange)) ; offset for first image matrix 
resuoff(O) = long(imagoff + total(simag)) ; offset for file result matrix 
for i = 0, npol - 2 do resuoff(i+l) = resuoff(i) + sresu(i) 

; Calculations 

for n = 0,npol-l do begin 
; create result matrix 
result = fltarr(100,nfra(n)+l) 

print, POLARIZATION CASE : ’+strtrim(string(n),l) 

; elapsed time for each frame 

time = leia_pte(filename, tempo’,n) 

etime = deltat(time) ; calculate elapsed time 

inclua, etime, result, 0 

; metdata 

info = leia_pte(filename, info’,n) 
time = leia_pte(filename, ’tempo ’,n) 

metoc = meteoro (info,time) ; interpolates metoc data 

; position data 

gpsship = metoc(0:2,*) ; GPS position of the ship 

camera = [1959.917,7034.517,27.439] ; lat. Ion, alt. of the camera in minutes 

auxones = make_array(l,nfra(n) Afloat, value=1.0) 

gpscam = camera#auxones 

position = Geography(gpsship, gpscam) 

; skindata 

info = leia_pte(filename, r mfo’,n) 
time = leia_pte(filename, ’tempo ’,n) 
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; interpolates skintemp data 



skin = skintemp(info,time) 

; include results in Result matrix 
inclua, metoc, result, 1 
inclua, gpscam, result, 10 
inclua, position, result, 13 
inclua, skin, result, 16 
; open file, localize pointer and save matrix 
escreva, result, filename, resuoff(n) 
endfor 
end 
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TAPE5SHIP.PRO 



.************** LXJG PANAGIOTIS KARAVAS HELLENIC NAVY***************** 
.********************** THESIS RF.SF ARCH* ************************ ****** 

;*** VALIDATION EVALUATION OF THE POLARIZATION VERSION OF SEARAD****** 
;PROGRAM NAME :Data.pro 

;Purpose : This program writes formatted input files tape5 files for the 

; specified geographical positions we choose 

pro tape5ship 

; define variables 

al=intarr(13) 

bl=fltarr(2) 

a2=intarr(6) 

b2=fltarr(5) 

a2a=fltarr(3) 

b2a=intarr(l) 

a2b=fltarr(3) 

a2c=intarr(3) 

b2c=strarr(18) 

a2cl=fltarr(6) 

b2cl=strarr(15) 

a2c2=fltarr(8) 

b2c2=fltarr(8) 

a2c3=fltarr(3) 

b2c3=intarr(5) 

a2d=intarr(4) 

a2dl=fltarr(l) 

b2dl=strarr(18) 

a2d2=fltarr(4) 

b2d2=fltarr(4) 

c2d2=fltarr(4) 

d2d2=fltarr(4) 

a3=fltarr(7) 

a3al=intarr(4) 

a3a2=fltarr(8) 

a3bl=intarr(l) 

a3b2=fltarr(6) 

b3=intarr(l) 

c3=strarr(l) 

a4=intarr(5) 

a5=intarr(l) 

;pick the file of interest 

basepath = ’c:\Thesis\data\basicV ; where is the B ASEN.PTE 

filename = pickfile(/read, path = basepath , filter = ^.pte’, title^elect file to Analyse) 
askl=0 

read, askl, prompt = OHOOSE POLARIZATION CASE HORIZONTAL = 0 VERTICAL =1 
UNPOLARIZED =2 : ’ 

if ((askl ne 0 )and (askl ne 1) and (askl ne 2) ) then return 
;get the information matrix for the specific image 
A=leia_pte(filename, r imagem’,askl) 

; dimensions & initial val 
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auxl = size(A) 

; Assign dimensions 
nima = 1 

if (auxl(O) eq 3) then nima = auxl (3) 

;choose desired frame to analyze 

print, THE NUMBER OF FRAMES FOR THIS POLARIZATION CASE IS:’ 
print, nima 

read, ask2, prompt = CHOOSE NUMBER OF FRAME TO ANALYZE’ 

;choose desired frame to analyze 
val=input(filename,ask 1 ,ask2) 

;get the name for the tape5 file 
N = ” 

read, N, prompt = 7PTE : Sequential number to store ship tape5 file (t .std) : ’ 

f=strtrim(string(fix(askl)),l) ;number that gives the polarization 

datapath = ’c:\Thesis\data\tape5Y+N+\shipV ;where the tape 5 files will be written 
;write the ship tape5 file 
;card 1 inputs 

al(0)=2 ; Mid-Latitude summer 
al(l)=2 ;slant path between to space 

al(2)=2 ; Radiance mode with solar/lunar scattered radiance included 

al(3)=l ; multiple scattering 

al(4)=0 ; default 

al(5)=0 ; default 

al(6)=0 ; default 

al(7)=0 ; default 

al(8)=0 ; default 

al(9)=0 ; default 

al(10)=0 ; default 

al(l 1)=0 ; normal operation(IM) 

al(12)=0 ; normal operation (NORT) 

bl(0)=273.15+val(5);Tbound 

bl(l)=0 ; SALB=0 assumes the surface is blackbody 

; card2 input 

a2(0)=3 ; Navy maritime extinction model 
a2(l)=0 ; Defaults to seasons of selected model’ 
a2(2)=0 ; Default to stratospheric background’ 
a2(3)=3 ; Air mass parameter 
a2(4)=0 ; No clouds or rain 
a2(5)=0 ; IVSA not used 
b2(0)=0 ; meteoroligical range 
b2(l)=val(3) ; current wind speed 

b2(2)=4.481 ; 24 hour average wind speed from metoc 2 data’ 
b2(3)=0 ; rain rate 
b2(4)=0 ; ground altitude 
;card3 input 

a3(0)=0.027 ; camera altitude 90 ft 27.4meters 
a3(l)=0.007 ; point altitude 

a3(2)=val(2) ; zenith angle of point 

a3(3)=0; range ,we do not need to specify it 
a3(4)=0 ; no need to specify it 

a3(5)=0 ; earth radius set to default 
b3=0 ; normal operation 
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a3(6)=val(4); wind azimuth relative to LOS 
c3=T’ ;seas witch 
;card3al input 

a3al(0)=l ; EP ARM defines card3a2 

a3al(l)= 2; MIE generated data base of aerosol phase functions for lowtran models 
a3al(2)=val(6);Julian day 
a3al(3)=0 ; sun as source of radiance 
;card3a2 input 

a3a2(0)=32.665 ; observer latitude (+ means North) 
a3a2(l)=l 17.242 ;observer longitude (0 to 360) 
a3a2(2)=0 
a3a2(3)=0 

a3a2(4)=val(7);Greenwich time in decimal hours 
a3a2(5)=val( 1 ) ; path azimuth 

if (val(7) gt 24.0) then begin 
a3a2(4)=a3a2(4)-24.0 
a3al(2)=a3al(2)+l 
endif 

;card4 input 

a4(0)=833 ; initial frequency (wavenumber cm-1 ) 
a4(l)=1250 ; final frequency (wavenumber cm-1 ) 
a4(2)=10 ; frequency increment (or step size) (cm-1) 
a4(3)=5 ; EFWHM 

;card5 input set to default mode i.e. end program 

; write input cards 

fname=datapa th+ ’s ’+N +f+ ’. std ’ 

openw,unit, fname, /get_lun 

printf, unit, format=’( H F",I4,12I5,F8.2,F7.2)’,al,bl 

printf, unit, format=’(6I5,5F10.3)’,a2,b2 

printf, unit, format=’(6F10.3,I5,F10.3,A5)’,a3(0:5),b3,a3(6),c3 

printf, unit,format=’(4I5)’,a3al 

printf ,unit,format= ’( 8F 1 0. 3) ’,a3 a2 

printf, unit, format=’(5I10)’,a4 

printf, unit, format=’(I5)’,a5 

free_lun,unit 

end 
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INPUT.PRO 



.************** LTJG PANAGIOTIS KARAVAS HELLENIC NAVY***************** 
.********************** THESIS RFSF. ARCH** ********************** ******* 

;***validation evaluation of the polarization version of searad****** 

;PROGRAM NAME :input.pro 
;PURPOSE: Provide input to SEARAD Tape5 files 
function input, filename, ask2,fn 
;read transmittance and path radiance 

shipout= c:\Thesis\data\out\’+fn+\ship\sVfn+2.std’;where the out file for the ship is 
tr=0.0 
pathr=0.0 

openr,unit,shipout/get_lun 

point_lun, unit, 1435 ;offset to get the zenith angle 
readf, unit,tr 

point_lun, unit, 1695 ;offset to get the zenith angle 

readf, unit,pathr 

free_lun,unit 

print, The path transmittance is :’,tr 
print, The path radiance is :’,pathr 

A=leia_pte(filename,Imagem’,2);get the unpolarized image 
pi = 3.1416 
;find hotspot 
hot=hotspot(A) 

;get the area of interest 
box=[0, 255,124,0] 

;get the horizon 
horizon=horizont(A) 

seam=element(A(*,*,ask2),hot,box,horizon) 

;convert digital level to thermal level and then to temperature 
dat=leia_pte(filename,’resultado’,2) 

Tl=dat(16,ask2+1) ;the temperature of the stack in C 

print, hot tmp’,Tl 

rad=radiance(T 1 ,tr) 

rad(0)=rad(0)+pathr 

print, The corected radiance is ’,rad 

Tref=invrad(rad) 

print, The corrected Temp is :’,Tref-273 
lrange=leia_pte(filename,lrange’,2) 
range=lrange( 1 ,ask2) 
level=lrange(0,ask2) 

B=isotherm(A(*,*,ask2),range, level) ;convert image matrix from Digital level to isotherm 

units 

tvscl,B,order=l 

h=isotherm(hot(0,ask2),range,level) ;convert hotspot digital level to IU 

print,hof,h 

print/max arr’,max(B) 

os=offset(h,Tref,2) ;Find the offset for the frame 
print, ’offset: ’,os 

reference point to aply geo data 
xref=hot( 1 ,ask2) 
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yref=hot(2,ask2) 
xdim=box(2)-box(0) 
ydim=box( 1 )-box(3) 
b=fix(xdim/9) 
c=fix(ydim/9) 

val=fltarr( 1 4,(b+ 1 )*(c+ 1)3) 
val(l l,*3)=os ;gives the offset 

val(12,*,2)=dat(15,ask2+l)*180/pi ;azimuth at the reference point 
val(13,*,2)=dat(14,ask2+l)*180/pi ;zenith angle at the reference point 

y=-9 

1=-1 

for i=0 ,c do begin 

y=y+9 ;steps of about 0.5 degrees for 7 deg lens 

x=-6 

for j= 0 ,b do begin 

x=x+9 ;steps of about 0.25 degrees for 7 deg lens 

;get the diference in azimouth and elevation 
dx=x-xref 
dy=y-yref 

daz i m=(-dx) *0 . 000977 4 
delev=dy*0.0004772 
azim=dazim+dat( 1 5 ,ask2+ 1 ) 
elev=delev+dat(14,ask2+l) ;zenith angle 

r=(27.439-7.3)/cos(180-(elev/pi)*180) ;compute the range from the zenith angle 
os=-0.0043*r + 12.162 

A(x,y)=temp(B(x,y),os, 2) ;matrix which contains temperatures 
1 = 1+1 

val(0,l,2)=A(x,y) * seam(x,y) temperature value calibrated 
;print,x,y 

val(l,l,2)=(elev/pi)*180 ;zenith angle at that point in deg 

val(2,l,2)=(azim/pi)* 1 80 ; true azimuth angle 

val(3,l,2)=dat(8,ask2+l)-val(2,l,2);wind direction relative to LOS 

val(4,l,2)=x 

val(5,l,2)=y 

endfor 

endfor 

vai(6,*,2)=dat(7,ask2+l) ;wind speed 

val(7,*,2)=dat(13,ask2+l) ;range to reference point in meteres 
val(8,*,2)=dat(9,ask2+ 1 ) ;sea surface temperatures 

info = leia_pte(fliename, r info’,2) 
time = leia_pte(filename, ’tempo'3) 
val(9,*3)=info(0)+90 

; transform time hhmmss.cc to ssssss.ss, make the corection + 8hours to make it GMT 
;and then divide by 3600 to take decimal hours GMT 

hour = (time(0,ask2)*3600.00+ 

time(l,ask2)*60.0+time(2,ask2)+time(3,ask2)/100.0+(8*3600.0))/3600.0 

val(10,*,2)=hour 

for i=0, 1 do begin 

A=leia_pte(filename/imagem’,i);get the polarized images 

dat=leia_pte(filename,resultado\i) 

lrange=leia_pte(filename,lrange’,i) 
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B=isotherm(A(*,*,ask2), range, level) ;convert image matrix from Digital level to isotherm 

units 

for j=0J do begin 

val(0,j,i)=B(val(4,j,2),val(5,j,2))*seam(val(4,j,2),val(5,j,2)) ;get the IU values of the polarized 
frames 

val( 1 :6,j,i)=val( 1 :6,j,2) 
val(8:ll,j,i)=val(8:l l,j,2) 
endfor 

val(7,*,i)=dat(13,ask2+l) ;range to reference point in meteres 
val(12,*,i)=dat(15,ask2+l)*180/pi;azimuth at the reference point 
val(13,*,i)=dat(14,ask2+l)*l 80/pi ;zenith angle at the reference point 
endfor 
return, val 
end 
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ISOTHERM.PRO 



.************** LTJG PANAGIOTIS KARAVAS HELLENIC NAVY******************* 
.sic************** ******* THESIS RF5?FAP PH* ************************ ************ 

;*** validation evaluation of the polarization version of searad****** 

;PROGRAM NAME risotherm.pro 
;PURPOSE: Convert Digital Level to Isotherm Units 
function isotherm, DL,range,level 
case range of 

5:begin 

a=0.00122I 

b=level-2.5 

end 

10:begin 

a=0.002442 

b-level-5 

end 

20:begin 

a=0.004884 

b=level-10 

end 

endcase 
IU=DL*a+b 
return, IU 
end 



INVRAD.PRO 

*************** LTJG PANAGIOTIS KARAVAS HELLENIC NAVY******************* 

*********************** thesis research** ************************ *********** 

;*** VALID ATI ON EVALUATION OF THE POLARIZATION VERSION OF SEARAD****** 
;PROGRAM NAME rinvrad.pro 

•JPURPOSE: Convert Radiance to Temperature, using Stefan Boltzmann’s Law 
;rad(0) radiance 

;rad(l) fraction of radiance in the 8-12 band 

function invrad,rad 

pi=3. 141592654 

s=5.66961 * 10.0 A (-8.0) 

e=l 

T=(pi*rad(0)/(rad(l)*e*s)) A 0.250 

retum,T 

end 
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RADIANCE.PRO 



LTJG PANAGIOTIS KARAVAS HELLENIC NAVY******************* 

. ****** ****** ********* sic THESIS RFSFARGH** *********************** ************ 
•*** VALIDATION EVALUATION OF THE POLARIZATION VERSION OF SEARAD****** 
;PROGRAM NAME iradiance.pro 

;PURPOSE: Calculate Radiance from Temperature, using Stefan Boltzmann Law. 

function radiance, T,tr 

pi=3. 141592654 

rad=fltarr(2,27) 

rad(0,0)=0.140 

for i=l,26 do begin 

rad(0,i)=rad(0,i- 1)+0.010 

endfor 

rad(l,*)=[0.0079053, 0.013023, 0.019962, 0.028858, 0.039754, 0.052613,$ 

0.067331,0.083750,0.10168,0.12091,0.14122,0.16239,0.18423,0.20653,$ 

0.22911,0.25183,0.27454,0.29712,0.31947,0.34150,0.36314,0.38432,0.40502,0.42518,$ 

0.44479,0.46382,0.48227] 

tl 1 = 8 . 0 * (T+27 3 . 0)/ 1 0000 

tl2= 12.0*(T +27 3 .0)/ 1 0000 
tl=[tll,tl2] 



res = INTERPOL(rad(l,*), rad(0,*), tl) interpolate. 

plot,rad(0,*),rad( 1 ,*),psym=2 

oplot,tl,res,psym=6 

fra=res(l)-res(0) 

print, Traction’, fra 

s=5 .6696 1 * 10.0 A (-8) 

e=l 

W=tr*fra*e*s*(T+273.0) A 4/pi 

;print,’radiance’,W 

result=fltarr(2) 

result(0)=W 

result(l)=fra 

return, result 

end 
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OFFSET.PRO 



*************** LTJG PANAGIOTIS KARAVAS HELLENIC NAVY******************* 
*********************** THESIS RESEARCH** ******************* **************** 

; * * * V ALID AT ION EVALUATION OF THE POLARIZATION VERSION OF SEARAD****** 
;PROGRAM NAME :Offset.pro 

;PURPOSE: Calculate Offset from input IU and Temperature 
function offset, TL, T, casep 
case casep of 
0 : begin 

A = double(1675) 

B = double(1348) 

C = double(l) 

end 

1 : begin 

A = double(1675) 

B = double(1348) 

C = double(l) 

end 

2 : begin 

A = double(6808) 

B =double(1511) 

C = double(l) 

end 

endcase 

os=TL - (A/(C*exp(B/T)- 1 )) 

return, os 

end 
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TEMP.PRO 



.************** LTJG PANAGIOTIS KARAVAS HELLENIC NAVY******************* 
*********************** THE STS RESEARCH************************************* 

;*** validation evaluation of the polarization version of searad****** 

;PROGRAM NAME :Temp.pro 

;PURPOSE: Calculate Temperature from input Offset and Temperature 

function temp, TL,os, casep 
auxl = size(TL) 
case casep of 
0 : begin 

A = double(1675) 

B = double(1348) 

C = doublet 1) 

;os=double(39) 

end 

1 : begin 

A = double(1675) 

B = doublet 1348) 

C = doublet 1) 

;os=double(39) 

end 

2 : begin 

A = double(6808) 

B = double(151 1) 

C = doublet 1) 

;os=double(12) 

end 

endcase 

w=alog(((A/(TL-os))+ 1 )/C) 

T=B/w-273 
return, T 
end 
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TAPE5.PRO 



.************** LTJG PANAGIOTIS KAlRAVAS HELLENIC NAVY******************* 
.********************** THESIS RESEARCH************************************* 
;*** VALIDATION EVALUATION OF THE POLARIZATION VERSION OF SEARAD****** 
;PROGRAM NAME :Tape5.pro 

;PURPOSE: This program writes formatted input files tape5 files for the 
; specified geographical positions we choose, for use by the polarized SEARAD 

; pro tape5 

; define variables 

al=intarr(13) 

bl=fltarr(2) 

a2=intarr(6) 

b2=fltarr(5) 

a2a=fltarr(3) 

b2a=intarr(l) 

a2b=fltarr(3) 

a2c=intarr(3) 

b2c=strarr(18) 

a2cl=fltarr(6) 

b2cl=strarr(15) 

a2c2=fltarr(8) 

b2c2=fltarr(8) 

a2c3=fltarr(3) 

b2c3=intarr(5) 

a2d=intarr(4) 

a2dl=fltarr(l) 

b2dl=strarr(18) 

a2d2=fltarr(4) 

b2d2=fltarr(4) 

c2d2=fltarr(4) 

d2d2=fltarr(4) 

a3=fltarr(7) 

a3al=intarr(4) 

a3a2=fltarr(8) 

a3bl=intarr(l) 

a3b2=fltarr(6) 

b3=intarr(l) 

c3=strarr(l) 

a4=intarr(5) 

a5=intarr(l) 

;pick the file of interest 

basepath = ’c:\Thesis\data\basicV ; where is the BASEN.PTE 

filename = pickfile(/read, path = basepath , filter = ^.pte’, title=$elect file to Analyse') 
A=leia_pte(filename, r imagem',2) 

; dimensions & initial val 
auxl = size(A) 



; Assign dimensions 
nima = 1 
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if (auxl(O) eq 3) then nima = auxl(3) 

;choose desired frame to analyze 

print, THE NUMBER OF FRAMES FOR THE UNPOLARIZED CASE ARE:’ 
print,nima 

read, ask2, prompt = CHOOSE NUMBER OF FRAME TO ANALYZE’ 
;choose desired frame to analyze 
;get the name for the tape5 file 
N = ” 

read, N, prompt = 7PTE : Sequential number to store data (t .std) : ’ 

datapath = ’c:\Thesis\data\tape5V+N+Y ;where the tape 5 files will be written 
val=inputcal(filename,ask2,N) 

j=0 

s=size(val(*,*,2)) 
for i=0,s(2)-l do begin 

if (val(0,i,2) GT 0.0) then begin 

j=j+l ;count the number of tape 5 files we will create 

;card 1 inputs 

al(0)=2 ; Mid-Latitude summer 
al(l)=3 ;slant path between to space 

al(2)=2 ; Radiance mode with solar/lunar scattered radiance included 

al(3)=l ; multiple scattering 

al(4)=0 ; default 

al(5)=0 ; default 

al(6)=0 ; default 

al(7)=0 ; default 

al(8)=0 ; default 

al(9)=0 ; default 

al(10)=0 ; default 

al(l 1)=0 ; normal operation(IM) 

al(12)=0 ; normal operation (NORT) 

b 1 (0)=273. 1 5+val(8,i,2) ;Tbound 

bl(l)=0 ; SALB=0 assumes the surface is blackbody 

; card2 input 

a2(0)=3 ; Navy maritime extinction model 
a2(l)=0 ; Defaults to seasons of selected model’ 
a2(2)=0 ; Default to stratospheric background’ 
a2(3)=3 ; Air mass parameter 
a2(4)=0 ; No clouds or rain 
a2(5)=0 ; IVSA not used 
b2(0)=0 ; meteoroligical range 
b2(l)=val(6,i,2) ; current wind speed 

b2(2)=4.481 ; 24 hour average wind speed from metoc 2 data’ 
b2(3)=0 ; rain rate 
b2(4)=0 ; ground altitude 

;card3 input 

a3(0)=0.027 ; camera altitude 90 ft 27.4meters 
a3(l)=0 ; point altitude 
a3(2)=val(l ,i,2) ; zenith angle of point 
a3(3)=0 ; range ,we do not need to specify it 
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a3(4)=0 ; no need to specify it 

a3(5)=0 ; earth radius set to default 
b3=0 ; normal operation 
a3(6)=val(3,i,2); wind azimuth relative to LOS 
c3=T’;seas witch 

;card3al input 

a3al(0)=l ; IP ARM defines card3a2 

a3al(l)= 2; MEE generated data base of aerosol phase functions for lowtran models 
a3al(2)=val(9,i,2);Julian day 
a3al(3)=0 ; sun as source of radiance 

;card3a2 input 

a3a2(0)=32.665 ; observer latitude (+ means North) 
a3a2(l)=l 17.242 ;observer longitude (0 to 360) 
a3a2(2)=0 
a3a2(3)=0 

a3a2(4)=val(10,i,2);Greenwich time in decimal hours 
a3a2(5)=val(2,i,2) ; path azimuth 

if (val(10,i,2) gt 24.0) then begin 
a3a2(4)=a3a2(4)-24.0 
a3a 1 (2)=a3a 1 (2)+ 1 
endif 

;card4 input 

a4(0)=833 ; initial frequency (wavenumber cm-1 ) 
a4(l)=1250 ; final frequency( wavenumber cm-1 ) 
a4(2)=10 ; frequency increment (or step size) (cm-1) 
a4(3)=5 ; IFWHM 

;card5 input set to default mode i.e. end program 
; write input cards 

g=strtrim(string(fix(j)),l) ; gives the number of the tape 5 file 
fname=datapath+ V+N+ l’+g+ \s td ’ 
openw,unit, fname, /get_lun 
printf, unit, format=’( M F M ,I4,12I5,F8.2,F7.2)’,al,bl 
printf, unit, format=’(6I5,5F10.3)’,a2,b2 
printf, unit, format= , (6F10.3,I5pl0.3,A5) , ,a3(0:5),b3,a3(6),c3 
printf, unit,format=’(4I5)’,a3al 
printf,unit,format=’(SF10.3)’,a3a2 
printf, unit, format=’(5I10)’,a4 
printf, unit, format=’(I5)\a5 
free_lun,unit 
endif 

endfor 

;get the name for the data file 

basename = DAT’+ strtrim(string(N),l)+’.PTE’ 



; create the DATA file 
vall=fltarr(6,3) 
value=fltarr(6,j,3) 
for k=0,2 do begin 
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vall(0: 
vall(2 
vall(3, 
vail (4, 
vail (5, 
m=-l 
for i=0. 



write 



l,k)=val(6:7,0,k);get the wind speed and distance to the ship 

k)=val(12,0,k) ;get the azimuth at the ship 

k)=val( 13,0,k) ;get the zenith at the ship 

k)=val( 1 l,0,k) ;get the offset for the temperature 

k)=j ;get the number of sea pixels that are analyzed 



s(2)-l do begin 

if (val(0 ,i,k) GT 0.0) then begin 
m=m+ 1 

vahae(0:5,m,k)=val(0:5,i,k) 

endif 



;choose useful elements of the val array to 



endfor 

endfor 

openw, unit, datapath+basename, /get_lun 

printf,unit,vall 

printf, unit, value 

close, unit 

free_lun, unit 

end 
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OUT.PRO 



.************** ltjg PANAGIOTIS KARAVAS HELLENIC NAVY******************* 

. ************* sic******** THESIS RFSF.aR PH** ********************* ************** 

Invalidation evaluation of the polarization version of searad****** 

;PROGRAM NAME :Out.pro 

;PURPOSE: This program reads formatted output files “out” from the polarized SEARAD, as well 
as experimental radiance, and provides an output file containing the polarization degrees and the 
radiance from the Experiment and SEARAD. 

pro out 



pi = 3.1416 
fn= ’ 

read, fn , prompt = ’CHOOSE FILE NUMBER : ' 
basepath=c:\Thesis\data\tape5Y+fn+Y 
basepath 1 =c:\Thesis\data\outY+fn+Y 
pn=2 

pns=strtrim(string(pn), 1 ) ;con vert pn to string 
vall=fltarr(6,3) 

filename=basepath+Daf+fn+\pte' 
openr, unit, filename, /get_lun 

readf, unit, vail ; wind speed ,range ,#of tape 5 files, offset,azimuth ,zenith at the ship 
nu=fix(val 1(5,0)) 
val=fltarr(6,nu,3) 
readf, unit, val 
u=fltarr(l l,nu) 
for i=0,nu-l do begin 

is=strtrim(string(i+l),l) ;convert i+1 to string 
file^asepathlVT’+fn+’J+is+'.std’ 
openr, unit,file,/getjun 

;offset to get the zenith angle 



; # polarization cases 



pointjun, unit, 980 
readf, unit,s 
u(0,i)=s 

pointjun, unit, 1020 
readf, unit,sl 
u(l,i)=sl 

pointjun, unit, 1432 
readf, unit,s2 
u(2,i)=s2 

pointjun, unit, 1810 
readf, unit,s3 
u(9,i)=s3 

pointjun, unit, 1876 
readf, unit,s4 
u(10,i)=s4 

pointjun, unit, 2900 
readf, unit,q 
u(3,i)=q 

pointjun, unit, 2958 



;offset to get the range 



;offset to get the path azimuth 



;offset to get subsolar lat 



;offset to get subsolar Ion 



;offset to get the unpolarized full range radiance 



;offset to get the full range unpolarized temperature 
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;offset to get the full range horizontal p radiance 



readf, unit,ql 
u(4,i)=ql 

point_lun, unit, 3872 
readf,unit,r 
u(5,i)=r 

point_lun, unit, 3958 ;offset to get the full range horizontal p temperature 

readf,unit,rl 

u(6,i)=rl 

point_lun,unit,3894 ;offset to get the full range vertical p radiance 

readf,unit,s 

u(7,i)=s 

point_lun,unit,3967 ;offset to get the full range vertical p temperature 

readf,unit,sl 

u(8,i)=sl 

free_lun,unit 

endfor 

uc=fltarr(6,nu) 

uc(0,*)=u(0,*) ;zenith angle 

uc(l,*)=u(3,*) ;radiance for the unpolarized case from searad 
uc(2,*)=searad(val(0,*,2));get the values of the experiment(unpolarized) 
camera = [1959.9,7034.52,27.439] ; lat. Ion, alt. of the camera in minutes 
sun=fltarr(3,nu) 

auxzeros = make_array(l,nu ,/float, value=0.0) 

sun=[u(9,*)*60,u(10,*)*60,auxzeros] ; lat. Ion, alt. of the sun in minutes 

auxones = make_array(l,nu,/float, value=1.0) 

gpscam = camera#auxones 

sunposition = Geography(sun,gpscam) 

sunazimuth=sunposition(2,*)* 1 80/pi 

uc(3,*)=u(2,*)-sunazimuth 

uc(4,*)=(val(0,*,0)-val(0,*,l))/(val(0,*,0)+val(0,*,l))*100;Degree of polarization from 
the experiment 

uc(5,*)=(u(5,*)-u(7,*))/(u(5,*)+u(7,*))*100;degree of polarization from searad 
;make a plot of the results 
a=[min(uc(2,*)),min(uc(l,*))] 
b=[max(uc(2,*)),max(uc( 1 ,*))] 
plot,uc(0,*),uc( 1 ,*)7ynoz,psym=2,$ 

titles ’experiment and searad unpolarized radiances ’,xtitle= Zenith angle’,ytitle= ’Radiance in 
W*m2/sr’,$ 

yrange=[min(a)-0.5,max(b)+2.5] 

oplot,uc(0,*),uc(2,*),psym=3 

valls=strtrim(string(vall),l) ;convert pn to string 

stat=fltarr(5) 

stat = moment(abs(uc(l,*)-uc(2,*)),sdev=sd) statistics ,mean,stdev of the temperature diff 

between searad and experiment 
stats=strtrim(string(stat),l) 
sds=strtrim(string(sd), 1 ) 

xyouts,min(uc(0, 5,c )),max(b)+3.4, , * searad temp.’ 
xyouts,min(uc(0,*)),max(b)+2.6,’. experiment temp.’ 
xyouts,(max(uc(0,*))+min(uc(0,*)))/2,max(b)+3. 4, ’offset: ’+valls(4) 
xyouts,(max(uc(0,*))+min(uc(0,*)))/2,max(b)+2.6,’range to ship: ’+valls(l) 
xyouts,(max(uc(0,*))+min(uc(0,*)))/2,max(b)+ 1 .8,’aver.abs. error: ’+stats(0) 
xyouts,(max(uc(0,*))+min(uc(0,*)))/2,max(b)+l,’stdv: ’+sds 
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window, 1 

statl = moment(abs(uc(4,*)-uc(5,*)),sdev=sdl) ;statistics ,mean,stdev of the temperature diff 

between searad and experiment 

stats 1 =strtrim(string(stat 1 ), 1 ) 

sds 1 =s trtrim( stri ng(sd 1 ) , 1 ) 

stats 1 ==strtrim(string(stat 1 ) , 1 ) 

sds 1 =strtrim(string(sd 1 ), 1 ) 

al=[min(uc(4,*)),min(uc(5,*))] 

bl=[max(uc(4,*)),max(uc(5,*))] 

plot,uc(0,*),uc(4,*),psym=3,$ 

title= ’experiment and searad degree of polarization’,xtitle= Zenith angle ’,ytitle= Percent 
Polarization ’,$ 

yrange=[min(al)-0.5,max(bl)+0.5] 
oplot,uc(0,*),uc(5,*),psym=2 
xyouts,min(uc(0,*)),max(bl)+0.5,’* searad % pol.’ 
xyouts,rmn(uc(0,*)),max(bl)+0.3,’. experiment % pol.’ 

xyouts,(max(uc(0,*))+min(uc(0,*)))/2,max(bl)+0.1,’aver.abs. error: ’+statsl(0) +W 

xyouts,(max(uc(0,*))+min(uc(0,*)))/2,max(bl)-0.1,’stdv: ’+sdsl 

statistics=fltarr(6) 

statistics(0)= val 1 s(4) ;offset 

statistics( 1 )=val 1 s( 1 ) ;range 

statistics(2)=stats(0) ;average absolute error 

statistics(3)=sds standard deviation 

statistics(4)=statsl(0) ;average absolute error in polarization 

statistics (5 )=sds standard deviation in polarization 

;write the results 

basename = ’OUT’+ fn+’PTE’ 

openw, unit, basepathl+basename, /get_lun 

printf,unit,statistics 

printf, unit, uc 

close, unit 

free_lun, unit 

end 
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APPENDIX F 



SEARAD INPUT AND OUTPUT FILES 

This Appendix contains three sample SEARAD files. The first one is a sample of 
a TAPE 5 file which is the input file for SEARAD, and provides all the necessary 
information for the calculation of the sea radiance. The second one is a sample of an OUT 
file generated by SEARAD which provides the values of unpolarized sea radiance, and its 
polarized components. Finally the third file is a sample OUT file for the case where we 
have a situation where we do not wish to calculate the sea radiation, and only wish to 
calculate the transmittance and the path radiance. This kind of OUT files was used in the 
calibration procedure. 
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SAMPLE TAPE5 (INPUT) FILE 



F232 1000000000 290.19 0.00 
3 0 0 3 0 0 0.000 5.814 4.481 0.000 0.000 

0.027 0.000 92.000 0.000 0.000 0.000 0 78.809 T 

1 2 99 0 

32.665 117.242 0.000 0.000 20.335 202.227 0.000 0.000 

833 1250 10 5 0 

0 

BRIEF EXPANATION OF THE TAPE 5 FILE 

Card 1 (Line 1): 2 = Model Midlatitude Summer, 3= Slant path to space, 2 = 
program execution in thermal radiance mode with solar/lunar scattered radiance included, 
1 = Program executed with multiple scattering, 290.12 = Boundary Temperature in K. 

Card 2 (Line 2): 3 = Navy Maritime extinction, 3 = Air mass character 3, 5.814 = 
Current wind speed (m/sec), 4.481= 24 hour average wind speed. 

Card 3 (Line 3): 0.0027 = Initial altitude in Km, 92 = Initial zenith angle (degrees) 
as measured from the initial altitude, 78.809 = Wind azimuth with respect to look angle, 
T = Sea Switch is True. 

Card 3A1 (Line 4): 1 = Method of specifying the solar geometry in card 3A2, 2 = 
MIE generated internal database of aerosol phase functions for the MODTRAN models, 
99 = Day of the year to specify the sun’s location in the sky. 
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Card 3A2 (Line 5): 32.665 = Observer latitude, 117.242 = Observer longitude, 
20.335 = Greenwich time in decimal hours, 202.227 = path azimuth. 

Card 4 (Line 6): 833 = Initial frequency in cm' 1 , 1250 = Final frequency in cm' 1 , 
10 = Frequency increment in cm' 1 , 5= Triangular Slit function. 
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SAMPLE OUT FILE 



***** SEARAD, A MODIFICATION OF LOWTRAN7 * * * * * 
DATE: 02/09/1999 TIME: 00:32:43.36 

THERMAL PLUS SOLAR RADIANCE MODE 
MULTIPLE SCATTERING USED 
MARINE AEROSOL MODEL USED 
WIND SPEED = 5.81 M/SEC 

WIND SPEED = 4.48 M/SEC, 24 HR AVERAGE 

RELATIVE HUMIDITY = 76. 1 1 PERCENT 
AIRMASS CHARACTER = 3.0 

VISIBILITY = .00 KM 

SLANT PATH TO SPACE 
HI = .027 KM 

HMIN = .000 KM 

ANGLE = 91.727 DEG 

FREQUENCY RANGE 

IV1 = 830CM-1 ( 12.05 MICROMETERS) 

IV2 = 1250 CM- 1 ( 8.00 MICROMETERS) 

IDV = 10 CM-1 

IFWHM = 5 CM-1 

EFELTER = 0 

SUMMARY OF THE GEOMETRY CALCULATION 



HI = 


.027 KM 


H2 = 


.000 KM 


ANGLE = 


91.727 DEG 


RANGE = 


.897 KM 


BETA = 


.008 DEG 


PHI = 


88.280 DEG 


HMIN = 


.000 KM 


BENDING 


= .001 DEG 
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LEN 



0 



SEA AT 290. 19 K REPLACES BLACK BODY BOUNDARY 
UPWIND = 83.345 DEG EAST OF LINE OF SIGHT 

SUMMARY OF OBSERVATION GEOMETRY 
BETA = .00807 DEG 

PATH AZIMUTH = 197.691 DEG EAST OF NORTH 

WIND AZIMUTH = 281 .036 DEG EAST OF NORTH 

RECEIVER LATITUDE = 32.665 NORTH OF EQUATOR 
RECEIVER LONGITUDE = 1 17.242 WEST OF GREENWICH 
FOOTPRINT LATITUDE = 32.657 NORTH OF EQUATOR 
FOOTPRINT LONGITUDE = 117.245 WEST OF GREENWICH 
SUBSOLAR LATITUDE = 7.057 DEG NORTH OF EQUATOR 

SUBSOLAR LONGITUDE = 124.444 DEG WEST OF GREENWICH 

VALUES SEEN FROM FOOTPRINT 

RECEIVER ZENITH ANGLE = 88.280 DEG 

RECEIVER AZIMUTH = 263.347 DEG WEST OF UP WIND 

SOLAR ZENITH ANGLE = 26.461 DEG 

SOLAR AZIMUTH = 84.829 DEG WEST OF UP WIND 

SOLAR SPECULAR TELT = 30.922 DEG ( 5.02 SIGMA, PROB = 3.284E-05) 

ZERO RANGE UNPOLARIZED VALUES 

SEA EMISSION = 19.85953 W M-2 SR-1 (AV. EMISS. .5997) 

SKY REFLECTION . = 10.26155 W M-2 SR-1 
SUN GLINT = .00001 W M-2 SR-1 

TOTAL RADIANCE = 30.12108 W M-2 SR-1 
BLACK BODY TEMP. = 1 1.9 C 

FULL RANGE UNPOLARIZED VALUES 
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SEA EMISSION 


= 15.11176 W M-2 SR-1 


SKY REFLECTION 


= 7.73348 W M-2 SR-1 


SUN GLINT 


.00000 W M-2 SR-1 



PATH TO FOOTPRINT = 8.44320 W M-2 SR- 1 ( A V. TRANS. .7515) 

TOTAL RADIANCE = 31.28845 W M-2 SR-1 
BLACK BODY TEMP. = 14.0 C 

HORIZONTAL VERTICAL (H-V)/(H+V) 

(W M-2 SR-1) (W M-2 SR-1) (%) 

ZERO RANGE POLARIZED VALUES 



SEA EMISSION 


8.80812 11.05140 -11.3 


SKY REFLECTION 


5.91030 4.35124 15.2 


SUN GLINT 


.00001 .00000 93.2 


TOTAL RADIANCE 


14.71843 15.40264 -2.3 



BLACK BODY TEMP. (C) -23.3 -21.3 

FULL RANGE POLARIZED VALUES 



SEA EMISSION 


6.70654 8.40522 -11.2 


SKY REFLECTION 


4.45081 3.28267 15.1 


SUN GLINT 


.00000 .00000 93.1 



PATH TO FOOTPRINT 4.22160 4.22160 .0 

TOTAL RADIANCE 15.37895 15.90950 -1.7 

BLACK BODY TEMP. (C) -21.3 -19.8 
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OUT FILE USED TO CALCULATE TRANSMITTANCE AND PATH 

RADIANCE 

***** SEARAD, A MODIFICATION OF LOWTRAN7 ***** 

DATE: 02/05/1999 TIME: 01:37:25.11 

THERMAL PLUS SOLAR RADIANCE MODE 

MULTIPLE SCATTERING USED 

MARINE AEROSOL MODEL USED 

WIND SPEED = 5.81 M/SEC 

WIND SPEED = 4.48 M/SEC, 24 HR AVERAGE 

RELATIVE HUMIDITY = 76.11 PERCENT 

AIRMASS CHARACTER = 3.0 

VISIBILITY = .00 KM 

SLANT PATH, HI TO H2 

HI = .027 KM 

H2 = .007 KM 

ANGLE = 90.825 DEG 

RANGE = .000 KM 

BETA = .000 DEG 

LEN = 0 

FREQUENCY RANGE 



IV1 = 
IV2 = 
EDV = 
1FWHM = 
EFILTER = 



830 CM- 1 ( 
1250 CM- 1 ( 
10CM-1 
5CM-1 
0 



12.05 MICROMETERS) 
8.00 MICROMETERS) 



SUMMARY OF THE GEOMETRY CALCULATION 



HI = .027 KM 

H2 = .007 KM 

ANGLE = 90.819 DEG 

RANGE = 1.400 KM 

BETA = .013 DEG 

PHI = 89.188 DEG 

HMIN = .007 KM 

BENDING = .000 DEG 

LEN = 0 

TBOUND SET TO . 10 K FOR MARINE SKY 
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INTEGRATED ABSORPTION = 140.66 CM-1 FROM 830 TO 1250 CM-1 

AVERAGE TRANSMITTANCE = .6651 

MAXIMUM RADIANCE = 4.957E-02 W M-2 SR-1 (CM-l)-l AT 830.0 CM-1 

MINIMUM RADIANCE = 1.813E-02 W M-2 SR-1 (CM-l)-l AT 1090.0 CM-1 

BOUNDARY TEMPERATURE = .10 K 

BOUNDARY EMISSIVITY = 1.000 

FILTERED RADIANCE = 1.150E+01 W M-2 SR-1 

BLACKBODY TEMPERATURE = -33.6 C 
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APPENDIX G 



RESULTS 

1. UNPOLARIZED RADIANCE PLOTS 

This Appendix contains a list of plots of sea radiance as a function of the camera 
zenith angle. Both the experimental, and the SEARAD calculated radiance, are plotted for 
various ship camera distances. 
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Figure Gl.l Experimental and SEARAD Sea Radiance as a Function of the Zenith Angle. 

Different Values of Radiance for the Same Zenith Angle Correspond to Different 
Azimuth Angles. The Ship Camera Slant Range is 1016 m. The Average Absolute Error 
is 2.96 Wm^sr' 1 and the Standard Deviation is 0.7 Wm^sr' 1 
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Figure G1.2 Experimental and SEARAD Sea Radiance as a Function of the Zenith Angle. 

Different Values of Radiance for the Same Zenith Angle Correspond to Different 
Azimuth Angles. The Ship Camera Slant Range is 1399 m. The Average Absolute Error 
is 2.07 Wm^sr" 1 and the Standard Deviation is 0.73 Wm^sr' 1 
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Figure G1.3 Experimental and SEARAD Sea Radiance as a Function of the Zenith Angle. 

Different Values of Radiance for the Same Zenith Angle Correspond to Different 
Azimuth Angles. The Ship Camera Slant Range is 1881 m. The Average Absolute Error 
is 2.98 Wm' 2 sr'' and the Standard Deviation is 0.55 Wm' 2 sr’' 
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Figure G1.4 Experimental and SEARAD Sea Radiance as a Function of the Zenith Angle. 

Different Values of Radiance for the Same Zenith Angle Correspond to Different 
Azimuth Angles. The Ship Camera Slant Range is 2201 m. The Average Absolute Error 
is 1.5 Wm^sr' 1 and the Standard Deviation is 0.63 Wm' 2 sr _1 
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Figure G1.5 Experimental and SEARAD Sea Radiance as a Function of the Zenith Angle. 

Different Values of Radiance for the Same Zenith Angle Correspond to Different 
Azimuth Angles. The Ship Camera Slant Range is 2704 m. The Average Absolute Error 
is 2.56 Wm' 2 sr 1 and the Standard Deviation is 0.73 Wm' 2 sr‘' 
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2. DEGREE OF POLARIZATION PLOTS 

This Appendix presents plots that compare, the polarized SEARAD degree of 



polarization, and the experimental degree of polarization. 



° Experimental Degree of Polarization (%) 
* SEARAD Degree of Polarization (%) 




Figure G2.1 Comparison of the Polarized SEARAD, and Experimental Degrees of 
Polarization. Different Values of Radiance for the Same Zenith Angle Correspond to 
Different Azimuth Angles. The Ship Camera Slant Range is 1016 m. The Average 
Absolute Error is 0.6 % and the Standard Deviation is 0.7 % 
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O Experimental Degree of Polarization (%) 
X SEARAD Degree of Polarization (%) 




Figure G2.2 Comparison of the Polarized SEARAD, and Experimental Degrees of 
Polarization. Different Values of Radiance for the Same Zenith Angle Correspond to 
Different Azimuth Angles. The Ship Camera Slant Range is 1399 m. The Average 
Absolute Error is 0.71 % and the Standard Deviation is 0.73 % 
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° Experimental Degree of Polarization (%) 
x SEAR AD Degree of Polarization (%) 




Zenith Angle (deg) 



Figure G2.3 Comparison of the Polarized SEARAD, and Experimental Degrees of 
Polarization. Different Values of Radiance for the Same Zenith Angle Correspond to 
Different Azimuth Angles. The Ship Camera Slant Range is 1881 m. The Average 
Absolute Error is 0.37 % and the Standard Deviation is 0.55 % 
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x SEARAD Degree of Polarization (%} 
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Figure G2.4 Comparison of the Polarized SEARAD, and Experimental Degrees of 
Polarization. Different Values of Radiance for the Same Zenith Angle Correspond to 
Different Azimuth Angles. The Ship Camera Slant Range is 2201 m. The Average 
Absolute Error is 0.72 % and the Standard Deviation is 0.63 % 
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° Experimental Degree of Polarization (%) 
* SEARAD Degree of Polarization (%) 




Figure G2.5 Comparison of the Polarized SEARAD, and Experimental Degrees of 
Polarization. Different Values of Radiance for the Same Zenith Angle Correspond to 
Different Azimuth Angles. The Ship Camera Slant Range is 2704 m. The Average 
Absolute Error is 0.48 % and the Standard Deviation is 0.73 % 
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